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ABSTRACT 


This  paper  discusses  several  approaches  to  scattering  from 
slightly  rough,  very  rough,  composite,  and  rough  spherical  surfaces. 

Average  incoherent  scattering  cross  sections  for  slightly  rough 
surfaces  are  obtained  using  a  perturbation  technique  for  a  perfect  con¬ 
ductor,  and  for  a  surface  material  which  is  homogeneous  (with  variable 
permeability,  |ir,  as  well  as  permittivity,  «r).  The  results,  correct 
to  the  first  order,  include  polarization  dependence,  and  particular 
solutions  with  curves  are  presented  for  the  vertical,  horizontal,  and 
circular  states.  They  are  compared  with  measured  data.  The  scatter¬ 
ing  process  for  this  class  of  surfaces  is  interpreted  physically. 

Three  optics  techniques  yielding  results  for  very  rough  surface 
scattering  cross  sections  are  presented,  and  it  is  shown  that  they  all 
give  the  same  solution  for  Gaussian  surfaces.  Both  homogeneous  and 
perfectly  conducting  surface  materials  are  treated,  and  polarization 
dependence  is  preserved.  The  restrictions  under  which  the  solutions 
are  valid  are  enumerated.  The  physical  interpretations  of  the  scatter¬ 
ing  mechanism  for  this  class  of  surfaces  are  discussed  and  shed  valuable 
insight  on  the  process.  Curves  are  shown  for  backscattering,  employing 
two  probability  models . 

Surfaces  consisting  of  both  very  rough  and  slightly  rough 
structures  together,  termed  composite  surfaces,  are  analyzed,  and  an 
heuristic  physical  derivation  of  their  scattering  cross  section  is  pre¬ 
sented.  Curves  for  backscattering  cross  section  are  shown,  and  com¬ 
parison  with  measured  results  confirms  the  validity  of  this  approach. 

Average  backscattering  cross  sections  for  rough  spherical  sur¬ 
faces  are  presented.  Both  the  coherent  and  incoherent  cross  sections 
are  given  for  slightly  and  very  rough  spheres.  These  results  have 
application  to  planetary  surface  radar  scatter  and  to  passive  satellite 
communications . 

The  large  volume  of  published  material  on  the  subject  is  reviewed. 
The  first  two  parts  dealing  with  slightly  rough  and  very  rough  surfaces, 
rather  than  retrace  detailed  mathematical  derivations  available  elsewhere, 
attempt  to  accomplish  five  goals:  (i)  to  compare  the  various  analyses 
already  available,  (ii)  to  elucidate  xplicitly  the  approximations  involved 
in  each  approach,  (iii)  to  explain  and  interpret  the  physical  process  be¬ 
hind  the  mathematics  producing  the  scattering,  (iv)  to  present  meaningful 
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closed  form  mathematical  solutions  and  curves,  many  of  which  have  not 
appeared  before,  for  the  scattering  cross  sections,  and  (v)  to  compare 
theory  with  measurements.  The  material  in  the  last  two  sections, 
treating  composite  and  rough  spherical  surfaces,  is  new  as  far  as  the 
Western  literature  is  concerned,  and  is  based  directly  on  the  results  of 
the  first  two  sections;  physical  interpretation  of  the  scattering  mechanism 
again  is  emphasized. 
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SCATTERING  FROM  SURFACES  WITH  DIFFERENT 
ROUGHNESS  SCALES;  ANALYSIS  /ND  INTERPRETATION 


I.  INTRODUCTION 

Serious  attempts  to  analyze  the  scattering  of  electromagnetic 
waves  by  rough  surfaces  and  interfaces  did  not  begin  until  about  1950. 
Research  on  this  subject  was  stimulated  to  a  large  extent  by  an  ac¬ 
cumulation  of  measured  data  on  radar  scattering  from  terrain  and  the 
sea  which  was  made  during  and  immediately  after  World  War  II. 

Another  stimulating  factor  was  the  significant  expansion  and  dis¬ 
semination  of  knowledge  in  the  area  of  applied  statistics  which  took 
place  in  the  decade  before  1950. 

Since  then,  analysis  of  rough  surface  scattering  has  made  use  of 
one  of  two  classes  of  surface  model.  The  first,  which  we  shall  term 
the  geometrical  model,  deals  with  surfaces  made  up  of  given  determin¬ 
istic  shapes,  but  arranged  or  distributed  in  a  random  fashion;  distri¬ 
butions  of  hemispherical  bosses  on  a  plane  is  an  example  of  such  a 
model  analyzed  by  Twersky  (1957).  The  second  which  we  shall  call 
the  statistical  model,  treats  the  roughness  height  itself  as  a  random 
variable  from  the  outset.  In  the  latter  case,  one  must  choose  or  specify 
the  surface  probability  distributions  or  correlation  coefficient.  The 
statistical  models  seem  more  physically  reasonable,  since  nature 
rarely  composes  a  surface  of  given  deterministic  shapes.  It  is  this 
second  class  of  rough  surface  models  which  concerns  us  here.  A 
review  of  the  literature  covering  many  of  the  rough  surface  scattering 
investigations  is  given  in  each  respective  section.  More  details,  and 
results  on  other  geometrical  and  statistical  models  can  be  found  in 
Radar  Cross  Section  Handbook  (1968). 

Statistically  rough  surfaces  can  be  divided  into  three  classes, 
distinguished  by  the  relative  roughness  height  compared  to  wavelength: 
slightly  rough  surfaces,  intermediate,  and  very  rough  surfaces.  The 
intermediate  class  has  roughness  features  which  lie  in  a  region  where 
neither  high  nor  low  frequency  approximations  apply.  Only  the  slightly 
rough  and  very  rough  classes  can  be  handled  mathematically.  Jn 
general,  different  techniques  must  be  employed  to  solve  these  two 
classes.  Much  information  can  be  obtained  about  the  scattering 
behavior  of  rough  surfaces  from  a  study  of  these  two  classes,  and 
it  is  not  unreasonable  to  expect  that  extrapolation  between  them  can 
give  valuable  insight  into  scatter  from  the  intermediate  class  of  rough 
surfaces. 
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Unfortunately,  as  with  all  new  fields  of  investigation,  errors  and 
oversights  have  occurred.  Soviet  and  Western  developments  in  this 
area  often  paralleled  each  other;  most  investigators'  articles  and 
results  went  unnoticed  by  their  foreign  counterparts.  Frequently, 
controversy  in  the  literature  over  mistakes  or  minor  details  tended 
to  obscure  the  positive  and  common  points  of  agreement.  Most 
serious  of  all,  the  emphasis  on  mathematical  detail  has  in  the 
majority  of  cases  tended  to  obscure  the  simple  physical  concepts 
and  processes  producing  the  scattering.  Also,  many  invetigators 
were  not  careful  in  stating  the  approximations  and  assumptions 
inherent  in  their  analyses;  this  has  often  led  to  misinterpretation 
and  misuse  by  others. 

We  shall  attempt  here  to  emphasize  the  physical  explanation 
of  the  scattering  mechanism  behind  the  mathematical  formulations 
and  solutions .  Care  is  taken  to  outline  the  approximations  and 
restrictions  which  apply  to  each  model.  Nonetheless  results  are 
presented  in  the  most  general  form  possible,  including  the  effects 
of  polarization  and  surface  constitutive  parameters.  Several  new 
results  and  explanations  are  offered  for  slightly  rough  surfaces, 
surfaces  consisting  of  composite  roughnesses  of  several  scales,  and 
rough  spherical  surfaces.  In  the  case  of  very  rough  surfaces  several 
existing  analyses  based  upon  different  physical  formulations  are 
compared  and  shown  to  give  identical  results. 

To  introduce  the  terminology  we  shall  use,  imagine  that  a 
perfectly  smooth  surface  becomes  slightly  rough.  The  average 
specularly  reflected  field  present  for  the  smooth  surface  decreases 
somewhat  as  the  roughness  grows;  in  addition  there  is  a  smaller 
amount  of  power  scattered  away  from  the  specular  direction  due  to 
the  presence  of  the  roughness.  The  former  we  shall  call  the  coherent 
component  while  the  latter  is  termed  the  incoherent  component  of  the 
scattered  field.  More  precisely,  if  one  surface  of  an  ensemble  is 
replaced  by  another  statistically  similar  member,  or  if  the  surface 
is  translating  slowly  with  respect  to  transmitter  and  receiver,  the 
coherent  component  is  distinguished  bv  the  fact  that  the  average  value 
of  the  complex  field  (existing  at  and  near  the  specular  direction)  is 
non-zero.  The  incoherent  field  fluctuates,  however;  its  phase  angle 
is  uniformly  distributed;  its  average  value  is  zero.  On  the  other  hand, 
the  average  scattered  incoherent  power  (or  incoherent  field  magnitude 
squared)  is  non- zero.  As  the  roughness  size  increases,  the  coherent 
component  decreases  and  the  incoherent  component  increases.  In  the 
very  rough  surface  limit,  all  of  the  scattered  power  is  incoherent 
although  it  may  be  sharply  beamed  in  the  specular  direction.  This  paper 
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is  concerned  mainly  with  the  incoherent  scattered  power,  except  for 
the  brief  discussion  below. 

The  coherent  field  scattered  from  a  slightly  rough  surface  is 
computed  in  the  same  fashion  as  it  is  for  a  smooth  surface.  In  con¬ 
trast  to  the  incoherent  component,  it  is  highly  dependent  upon  the 
shape  of  the  illuminated  area.  It  exhibits  the  familiar  lobed  structure 
of  a  flat  plate,  for  example,  while  the  incoherent  component  does  not. 
Being  concentrated  around  the  specular  direction,  it  can  often  be  ignored 
in  other  directions.  Experimentally,  one  can  isolate  the  component  by 
averaging  first  and  then  squaring  the  received  voltage.  The  specular 
reflection  coefficients  for  the  slightly  rough  surface  may  be  written  in 
terms  of  the  Fresnel  coefficients  for  a  smooth  surface  and  a  factor 
which  accounts  for  roughness. 

This  modification  is  discussed  in  Radar  Cross  Section  Handbook 
(1968)  and  can  be  determined  for  a  Gaussian  surface  height.  The  result¬ 
ing  electric  field  reflection  coefficients  for  that  case  are 


(1) 


R 


R 11  (®i)  -  R11  (^i)  e 


■2ko  hZcosZ0j 


R?(6i)  =  Rx(0i)e 


-2ko  h2  cos2  0j 


where  RM  (0^  and  R^Gj)  are  the  usual  Fresnel  reflection  coefficients  for 
the  smooth  surface  and  are  given  in  Eq.  (35).  The  quantity  0j  is  the 
angle  of  incidence  from  normal,  kD  =  2tt/X.  is  the  wave  number,  and  h 
is  the  rms  roughness  height.  As  can  be  seen  from  the  above,  the  effect 
of  roughness  on  the  coherent  scattered  field  is  contained  in  the  ex¬ 
ponential  factor.  It  is  this  factor  which  quantitatively  expresses  the 
well-known  Rayleigh  criterion.  Clearly,  for  kQh  <  l/4,  the  coherent 
component  will  be  not  significantly  different  from  that  of  a  smooth  surface, 
and  we  call  the  surface  "slightly  rough".  On  the  other  hand  if  kQh  cos  0j>  1, 
the  coherent  component  will  be  insignificant,  and  we  call  the  surface 
"very  rough".  The  coherent  component  disappears  as  the  roughness 
height,  h,  increases. 

II.  SLIGHTLY  ROUGH  PLANAR  SURFACE 
A.  Introduction 

In  this  section  we  will  discuss  an  important  class  of  random  rough 
surfaces:  those  having  a  small  scale  of  roughness  whose  rms  roughness 
height,  h,  is  much  less  than  wavelength.  In  particular  we  shall  obtain 
expressions  for  the  incoherent  average  scattering  cross  sections  per 
unit  surface  area,  <r°,  for  various  polarization  combinations.  As 
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mentioned  in  the  introduction,  this  cross  section  accounts  for  most  of 
the  scattered  energy  everywhere  except  in  the  neighborhood  of  the 
specular  direction.  At  and  near  this  direction,  the  coherent  component 
is  strong  and  cannot  be  neglected. 

Slightly  rough  surface  models  have  been  employed  in  the  past 
based  upon  the  tangent  plane,  or  Kirchhoff  approximation.  The  first 
such  analysis  was  that  of  Davies  (1954).  It  has  been  discussed  also  by 
other  workers,  and  a  readable  treatment  is  found  in  Beckmann  and 
Spizzichino  (1963  Section  5.3).  It  should  be  mentioned,  however,  that 
those  models  are  often  inapplicable  to  natural  surfaces.  The  reason  is 
that  the  tangent  plane  approximation  restricts  one  to  surfaces  having 
2irpA,  >  >  1,  where  p  is  the  local  radius  of  curvature  of  the  surface 
at  all  points.  This  requirement,  along  with  the  defining  requirement 
2trh/\  <  <  1  generally  describes  a  very  smooth,  undulating  surface  with 
hills  much  less  than  wavelength  in  height  and  hundreds  of  wavelengths 
apart;  the  maximum  slopes  of  these  hills  are  infinitesimal.  What  is 
more,  the  model  fails  in  the  low  frequency  limit,  since  wavelength 
eventually  becomes  larger  than  surface  radii  of  curvature.  Also,  the 
model  shows  no  polarization  dependence  for  backscatter,  as  is  typical 
of  models  based  upon  the  tangent  plane  approximation.  Experimental 
evidence  shows  that  radar  backscatter  does  depend  upon  polarization; 
this  will  be  discussed  again  later. 

For  the  above  reasons,  we  shall  avoid  use  of  the  Kirchhoff 
approximation  for  the  slightly  rough  surface  model.  Instead,  we  will 
resort  to  a  much  less  restrictive  model,  but  interestingly  enough,  a 
model  originally  formulated  before  that  of  Davies  by  Rice  (1951). 

Rice's  formulation  was  later  developed  by  Peake  (1959)  who  derived 
the  scattering  cross-sections.  This  model  is  based  on  a  perturbation 
approach.  Since  the  tangent  plane  assumption  is  not  employed,  all 
of  the  accompanying  approximations  such  as  neglect  of  multiple  scat¬ 
tering  and  shadowing  are  avoided.  The  model  is  therefore  valid  for 
grazing  incidence  and  scattering  angles,  and  exhibits  meaningful 
polarization  dependence.  Most  important,  it  is  valid  in  the  low  fre¬ 
quency  limit  and  can  apply  to  any  of  the  class  of  very  rough  surfaces 
to  be  examined  later  when  the  frequency  is  low  enough. 

Valenzuela  (1967)  recently  extended  the  same  technique  to  obtain 
expressions  for  the  second  order  perturbation  correction  terms  for 
backscatter  from  dielectric  and  perfectly  conducting  surface.  Only  the 
first  order  terms  of  the  solution  are  obtained  here. 

The  following  are  the  restrictions  on  the  class  of  surfaces  to  which 
the  model  applies: 
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(a)  k0£,(x,  y)  <  1.0,  i.  e. ,  roughness  height  is  small  compared 
to  wavelength,  \;  £,(x,  y)  is  the  height  of  the  rough  surface  above  the 
x-y  plane,  which  is  taken  to  be  the  mean  surface  plane. 

(b)  8t,/3x,  a^/ay  <  1.0,  i.  e.  ,  surface  slopes  are  relatively 

small. 

(c)  <(di/dx2)>  =<(di/dyi)>,  i.  e.,  the  roughness  is  isotropic. 

Here,  the  brackets  <  >  indicate  average  over  an  ensemble  of  surfaces. 
This  restriction  is  not  essential  to  the  solution,  but  is  employed  here 
for  simplicity. 

(d)  L  >>  f,\,  i.  e. ,  the  dimension,  L,  of  the  illuminated  area 
is  much  greater  than  either  the  roughness  correlation  length,  l,  or 
the  wavelength. 

Since  the  tangent  plane  approximation  is  not  employed,  the  mutual 
interaction  effects  of  shadowing  and  multiple  scattering  are  not  over¬ 
looked.  Hence  the  slightly  rough  surface  solution  includes  these 
mechanisms.  In  fact,  mutual  interaction  between  neighboring  hills 
is  an  important  contributor  when  the  distance  between  them  is  less 
than  a  wavelength. 

We  shall  obtain  first  results  for  the  average  incoherent  scat¬ 
tering  cross  section  per  unit  surface  area  for  a  surface  of  homogeneous 
material  having  relative  permeability,  pr,  and  permittivity,  £r.  In  this 
respect  our  results  differ  from  the  analysis  of  Rice  (1951),  Peake  (1959), 
and  Valenzuela  (1967) ,  who  treat  only  dielectric  surfaces.  The 
quantities  pr  and  €r  may  be  real  or  complex,  representing  a  lossless 
or  lossy  material.  Their  magnitudes  may  be  greater,  equal  to,  or 
less  than  unity,  so  that  the  important  class  of  plasma  media  can  be 
included.  As  a  special  case,  the  imaginary  part  of  €r  can  approach 
infinity,  yielding  the  proper  results  for  a  perfectly  conducting  surface. 

The  notation  of  Rice  (1951)  ,  Peake  (1959) ,  and  Valenzuela  (1967) ,  will 


Also  implied  in  the  definition  of  scattering  cross  section  is  the  restriction 
that  the  observation  point  must  be  in  the  far  field  of  the  surface.  In  the 
absence  of  roughness,  this  requirement  means  that  the  distance  to  the 
observation  point,  R,  must  satisfy  the  criterion  R  >  2L2/\.  However, 
it  can  be  shown  (see  Barrick,  1965)  that  when  one  is  considering  the 
incoherent  scattered  power  from  a  roughened  surface,  this  requirement 
reduces  to  R  >  2f2A*  where  i  is  the  roughness  correlation  length. 

Physically  this  means  that  one  can  be  considerably  closer  to  a  rough 
surface  than  a  smooth  one,  since  f  <<  L  by  (d)  above. 
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be  employed  as  much  as  possible  here,  so  that  readers  wishing  a  more 
detailed  derivation  can  refer  to  that  work.  The  solution  for  a  horizintally 
polarized  incident  plane  wave  will  be  obtained  here.  Similar  results 
for  vertical  incident  polarization  will  be  then  written  from  inspection 
by  interchanging  the  roles  of  E  and  H,  «r  and  |ir.  From  these  results, 
the  pertinent  values  for  the  circular  polarization  states  and  for  linear 
polarization  states  of  arbitrary  orientation  will  be  given. 

B.  Analysis 

The  time  convention  used  here  will  be  e”lw*.  The  analyses  of 
Rice,  Peake,  and  Valenzuela  are  compatible  with  this  if  one  replaces 
i  there  by  -i,  since  they  use  the  positive  time  convention.  Let  us  first 
consider  one  particular  rough  surface  from  an  ensemble  of  statistically 
similar  rough  surfaces.  For  convenience,  assume  the  surface  to  be 
square,  of  length  L,  centered  about  the  origin  and  with  its  mean  surface 
coinciding  with  the  x-y  plane  (Fig.  1).  Later  in  the  analysis,  the 
surface  length  will  be  permitted  to  approach  infinity;  the  square  shape 
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of  the  surface  is  irrelevant  to  the  results,  but  is  a  useful  artifice  in 
performing  the  derivation.  Under  these  conditions,  the  surface  height 
can  be  expanded  in  the  following  Fourier  series: 

CO 

(2)  U*.  y)  =  £  £  P(m,  n)  eia<™+“y>  , 

m, n=-» 

,  _  2it 

where  a  =  —  . 

L 

The  real  nature  of  the  surface  height  requires  that  P(-m,  -n)  = 

P*  (m,  n)  where  the  asterisk  denotes  complex  conjugate.  The  trick  in 
simplifying  the  analysis  is  to  expand  the  scattered,  or  perturbed,  field 
into  a  Fourier  series  with  eigenfunctions  e*a(mx+ny).  Physically  this 
means  representing  these  fields  as  a  superposition  of  plane  waves  in 
all  directions.  Remembering  that  the  incident  field  is  horizontally 
polarized  in  the  y-direction  as  shown  in  Fig.  1,  the  total  E -field  in 
free  space  above  the  surface  is  then  written  as  follows: 


(3a)  Ej  =  Eax 

(3b)  Ey  =  E1  +  Er  +  E®, 

(3c)  Ez  =  Ez  * 


where  El  is  the  incident  field,  and  Er  is  the  reflected  field  from  a 
perfectly  smooth  homogeneous  surface.  The  terms  with  the  super¬ 
script,  s,  are  then  due  solely  to  the  presence  of  roughness,  and  may 
be  considered  the  incoherent  scattered  or  perturbed  fields.  These 
quantities  are  given  as 


(4a) 


(4b) 


E1  =  Eoe-<r«x-hrVz  Er  =  R  £o  e-crax-crVz 

y  y 


00 


E®  =  yy  Amneia(mx+nV>  •  eib<m-'1>z. 


m, n= ■» 
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r 


(4c) 


(4d) 


CO 


_  B^^eia(mx+ny)  #  eib(m,  n)  z  ^ 


,s 

JY  /  .  /  .  ~mn 

m.  n  =  -® 


CO 


C  ia(mx+ny) 
'-mn  e 


>ib(mJ  n)  z 


m,  n=-oo 


In  these  equations,  cr  =  -ikQ  =  iu >J  cQ\iQt  a  =  sin  0j,  Y  =  cos  0j, 
and  E°  is  the  electric  field  strength  of  the  incident  field.  R  is  the 
Fresnel  reflection  coefficient  for  horizontal  polarization,  and  is 
given  by 


(4e) 


R 


1  - 


1  + 


T  V' 

W-r°~  v 


Hr°-Y 


where  T  =  ik and  Y' 


1 


The  function  b(m,  n)  is  not  independent  of  the  wavenumbers  in  the 
x  and  y  direction,  and  is  dictated  by  the  requirement  that  each  term  of 
the  above  series  must  satisfy  the  wave  equation.  This  results  in  the 
following  relationship: 


(5) 


b(m,  n) 


a  (mz  + 


n2) 


4 


+  a  (mz  + 


n  ). 


The  objective  is  to  derive  Amn,  Bmn,  and  Cmn  in  terms  of 
P(m,  n)  and  all  of  the  above  parameters.  This  is  done  by  a  standard 
perturbation  technique  and  the  details  are  similar  to  those  of  Rice, 
Peake,  and  Valenzuela.  Briefly,  these  three  coefficients  are  obtained 
by  expanding  all  quantities  possible  into  series,  each  term  of  which 
involves  a  higher  power  of  the  smallness  parameters,  which  are  taken 
to  be  kQ4,  9t,/8x,  and  9t,/9y.  By  assumptions  (a)  and  (b) ,  these 
quantities  are  small,  and  consequently  only  the  first  term  of  these 
series  need  be  retained.  The  above  three  unknown  coefficients  (as 
well  as  three  other  unknown  coefficients,  Gmn,  Hmn,  Imn  which 
represent  the  scattered  fields  inside  the  surface  material)  are 
obtained  from  six  equations,  four  of  which  come  from  the  boundary 
conditions  for  the  E  and  H  fields  at  the  interface  z  =  £,,  and  two  of 
which  arise  from  the  divergence  conditions  on  each  side  of  the  interface. 
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The  three  desired  coefficients  Amn,  Bmn,  and  Cmn  are  then 
correct  to  the  first  order,  i,  e.  ,  the  neglected  terms  are  at  least  of 
order  kQC,  9&/9x,  or  94/9y  with  respect  to  the  result  shown.  The 
terms  of  the  next  higher  order  are  derived  by  Rice  (1951)  and 
Valenzuela  (1967)  for  perfectly  conducting  and  dielectric  surfaces, 
but  they  are  not  given  here.  Also  it  should  be  noted  that  each  of 
these  coefficients  is  directly  proportional  to  P(m-v,n),  the  m-v, 
n-th  coefficient  in  the  Fourier  expansion  for  the  surface  itself. 

( v  5  kQ  sin  9j/a,  and  is  taken  to  be  an  integer  here,  i.  e.  ,  the 
possible  directions  of  incidence  are  quantized  for  the  moment. 

However,  in  the  limit  as  a  =  2it/L  -*  0  with  L  —  ®  ,  v  becomes 
large,  and  thus  it  is  possible  to  approximate  it  by  the  nearest  integer 
with  negligible  error.  ) 

Up  to  now,  no  restrictions  have  been  placed  on  the  nature  of 
the  surface.  The  surface  may  in  fact  be  deterministic  or  periodic. 

For  example  in  the  case  of  a  one -dimensional  sinusoidal  surface  of 
period  L,  the  expansion  for  the  surface,  £(x),  degenerates  to  the 
situation  wher<*  P(j,  k)  3  0  for  j  ^  ^  1  and  k  ^  0;  this  importam 
case  was  solved  long  ago  by  Lord  Rayleigh  (Dover,  1945).  To  gain 
insight  into  the  scattering  process  for  the  slightly  rough  surface,  it 
is  instructive  to  point  out  a  fact  about  this  surface.  The  vanishing 
of  all  but  two  of  the  coefficients  P(m-v,  n)  means  that  the  coefficients 
for  the  scattered  field  expansions  of  Eqs.  (3b),  (3e),  i.  e. ,  Amn, 

Bmn»  an<*  Cmni  all  vanish  except  when  n  =  0  and  m  =  v  *  1.  In  other 
words,  the  scattered  field  consists  of  only  two  plane  waves,  propa¬ 
gating  in  two  directions  determined  by  the  direction  of  incidence  and 
the  wavelength-to-surface  period.  These  waves  are  both  in  the  plane 
of  incidence  (i,  e. ,  x-z  plane)  ,  and  the  scattering  angle  is  determined 
by  the  equation  sin  9g  =  sin  9^  ^  X/L.  It  is  important  to  note  that  these 
waves  exist  physically  only  when  |  sin  0S  |  =  |  sin  -  X./l|  <  1;  when 
k/L  is  such  that  this  inequality  is  not  satisfied,  the  propagation  number, 
b(m,  n)  is  pure  imaginary,  indicating  an  evanescent  or  non-radiating 
mode.  For  period  much  larger  than  wavelength,  scattering  takes  place 
close  to  the  specular  direction.  As  \/L  becomes  larger,  the  two 
scattered  plane  waves  diverge  from  the  specular  direction  until  they 
reach  grazing  (|  sin  9g  |  =  1),  at  which  point  the  plane  waves  are  cut 
off.  Hence  the  scattering  directions  are  determined  not  by  the  ampli¬ 
tude  of  the  sinusoidal  surface  (so  long  as  it  satisfies  the  "slightly  rough" 
criterion)  but  by  its  period  compared  to  wavelength.  This  behavior  is 
also  exhibited  by  the  random  rough  surface  considered  below. 

Before  discussing  the  scattering,  let  us  discuss  the  parameters 
involved  in  the  statistical  description  of  the  surface,  £,(x,  y).  Davenport 
and  Root  (1958)  employ  Rice’s  earlier  work  (1944,  1945)  to  show  that 
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a  random  but  periodic  function  of  period  L  can  be  expanded  in  a 
Fourier  series  such  as  Eq.  (2),  where  the  coefficients  P(m,  n)  are 
random  variables.  Going  further,  they  show  that  these  coefficients 
become  uncorrelated  in  the  limit  as  L  -*  oo  .  Practically,  it  is 
necessary  only  that  L  >  >  I ,  i.e.,  the  total  surface  length  must  be 
much  greater  than  the  roughness  correlation  length;  this  is  exactly 
restriction  (d)  above.  Now,  recalling  that  the  mean  value  of  the  sur¬ 
face  is  zero,  these  conditions  become. 


(6a)  <P(m,  n)>  =  0 


(6b) 


<P(m,  n)  P(u,  v)>  *  J 


0  for  u,  v  =  m,  n 


W(p,q)  for  u,  v  =  m,  n 


where  p  *  am  =  2irm/L  and  q  =  an  =  2irn/L.  The  function  W(p,  q)  is 
physically  the  roughness  spectral  density  of  the  surface,  and  p,  q  are 
radian  wave  numbers.  Using  these  equations,  the  following  relation¬ 
ships  are  obtained. 


(7a) 


<l‘  (x.  y)>  =  £  <P(m,n)  P(u.v)>  eiax(m-u)  +  iay‘n-v) 


m,n,u,v 


=  ^  <P(m,n)  P( -m,  -n)>  —  J  dm  ^  dn  jj  W(p,q) 


-00  -00 


/  m,n 


■irr 


W(p,  q)  dpdq  =  h2  , 


-00  -00 


(7b) 


<C(x,y)Ux',y't>  =  V  <P(m,n)P(u,v)>.iamx  +  iaUX'+iany+ iavy' 


m,n,u,v 


=  7  <P(m,n)P(-m,-n)>eiam(x-x')  +  ianly-y') 

J 

m,n 

— ’  ^  ^  W(P*  q)  elpTx  q  ydpdq  =  h2R(TX,  Ty) 

v  -00  -00 
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where  tx  =  x-x1,  and  Ty  =  y-y*.  The  quantity  h2 ,  as  mentioned  earlier, 
is  the  surface  mean  square  height,  and  R(tx,  Ty)  is  the  surface  height 
correlation  coefficient.  Relationship  (7b)  merely  states  that  the  rough¬ 
ness  height  spectral  density  and  the  surface  height  correlation  function 
are  Fourier  transforms. 

Now,  let  us  determine  the  horizontal  (h)  and  vertical  (v)  com¬ 
ponents  of  the  plane  wave  scattered  field  in  a  given  direction  represented 
by  the  angles  0g,  4>g;  first  it  is  necessary  to  relate  the  wave  numbers 
am  and  an  of  the  exponentials  of  Eqs .  (4b)-(4d)  to  the  propagation  con¬ 
stant  of  the  plane  waves  in  spherical  coordinates.  This  relationship 
is  obvious: 

(8)  am  =  kQ  sin  0S  cos  4>s»  and  an  =  kQ  sin  0S  sin  4>s. 


Then  the  horizontal  and  vertical  components  of  the  scattered  plane 
wave  in  the  direction  corresponding  to  m,  n  are 


(9) 

and 

i%) 


Ej(m.n)  =  (-Ann  .in  +  BmnCOS  +.)  .ia‘mx  +  T>  +  lb(m’ 

E“(m,n)  =  (Amncos  0S  cos  <Ca  +  Bm„cos  6S  sin  <|>s 
-  Cm„  sin  0.)  eia‘mx  +  “*>  +  ib(m>  n)z. 


Let  us  redefine  these  relationships  in  terms  of  quantities 
(ahh  an<*  avh)  w^ich  are  directly  proportional  to  scattering  matrix 
elements,  namely 


(10.)  Ej(m.n)  =  2k0  E°  ahh  cos  0(  P(m-v,n)  eiatmx  +  ny>  +  ib*m-  n>z. 
and 

( 10b)  E*(m,  n)  =  2k0  Ey  ovh  cos  0,  P(m-v.n)  cia(™*  +  ny>  +  ib(m,n)*. 


The  quantities  and  0y ^  are  determinable  from  Eqs.  (9)  and  the 
expressions  for  Amn»  Bmn»  Cmn.  They  will  be  given  in  simplified  form 
subsequently,  along  with  expressions  for  the  other  two  elements  a^y  and 
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avv  for  vertical  incident  polarization.  Note  that  the  random  nature  of 
the  rough  surface  is  reflected  explicitly  in  the  surface  coefficients 
above,  P(m-v,  n) .  The  factors  and  avh  are  nonrandom,  known 

functions.  At  this  point,  the  average  scattering  cross  sections  can  be 
determined.  First,  let  us  find  the  increment  in  average  field  intensity 
in  all  of  the  plane  waves  between  mi  and  m2,  ni  and  n2.  For  arbitrary 
polarizations,  y,  5,  this  is  given  by 


(11) 


m2  na 

A<|e®  (m,n)|2>  =  4k02  lEg^cos2^  ^  ^ 

mi  ni 


<P(m-v,  n)  P*(m-v,n)> 


where  use  is  made  of  the  orthogonal  properties  of  the  coefficients,  as 
stated  in  Eq.  (6b).  If  the  length  of  the  surface,  L,  is  large  compared 
to  wavelength,  X  ,  there  are  many  such  plane  waves  between  mi  and 
m2  (as  well  as  between  ni  and  n2  ).  The  terms  in  the  summation  are 
then  nearly  constant  for  mi  <  m  <  m2  (as  well  as  ni  <  n  <  n2  )  for  m2 
sufficiently  near  mi,  and  the  summation  can  be  rewritten 


(12)  A<  |E®(m,n)|2>  =  4k02  |Eg  |2  cos20i|aY&|2 

<P(m-v,  n)  P*(m-v,  n)>  (m2-mi)  (n2-ni) 


2  w  I  TT.O  |  2  ...2  q  l~  12 


=  4k0  jJ  |E&  l  2  cos  I  ayB  i 

W(am-  a  v>  an)AmAn, 


where  Am  =  m2  -  mi  and  An  =  n2  -  ni  . 

Using  the  above  expression,  the  increment  of  intensity  per  unit 
of  solid  angle,  fis,  is  expressed  in  terms  of  the  Jacobian,  AmAn/Afis  • 


(13) 


A  ^  Ey  ( bg  »  $  s )  2  ^ 
A  fig 


4k0TT2 

L2 


E°  | 2  cos2 


2 


W[  kQ(sin  6S  cos  4>s  " 


sin  Bj),  kQ  sin  bg  sin  <t>sJ 


Arr 

~Af 
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Noting  that  ft  s  =  sin  6S  dOs  dcj>s,  the  Jacobian  is 


(14) 


AmA  n 
Aft  s 


cos  0S 


Hence  the  average  field  intensity  at  the  receiver,  or  observation  point, 
is  the  quantity  in  Eq.  (13)  multiplied  by  the  solid  angle  subtended  by  the 
surface  from  the  observation  point.  This  is 


Aft  <5  = 


cos  0S 

1 - 


R 


so  that 


,i5)  <|E»|S.a<I*yw,,*'>I‘>  lV°^ 

»  Aft  s  R 

The  average  scattering  cross  section  per  unit  surface  area  at  the 
observation  point  and  for  incident*  scattered  polarization  states  5,  y  is 
given  by 


or 


(16) 


^y 6  (ys»  ^s) 


<  I  E§  1 2>  4ttRz 

l2  i=sTa 


°y6(es.<t>s)  =  4uko4  cos2  Oi  cos2  ys  la  y5  |2 

W[  kQ(sin  0S  cos  4>s  -  cos  0j),  kQ  sin  0S  sin  4>s]  . 


Restriction  (c)  assumed  that  the  surface  roughness  is  isotropic. 
This  means  that  the  height  correlation  function  R(tx,  Ty)  is  a  function 
only  of  the  separation,  r,  between  the  surface  points  x,  y  and  x\y', 

(i.e.,  r  =  J(x-x')2  +  (y-y1)*  =  Jtx2  +  t  ^2  )  •  Hence,  define 
(17a)  p(r)  =  r<tx>  Ty), 


applicable  when  the  roughness  is  isotropic.  This  means  that  the  rough¬ 
ness  spectral  density,  W(p,  q),  is  also  symmetrical  in  p  and  q.  Define 
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t  =  Jp2  +  q2  and 
(17b)  w(t)  =  W(p,q). 


Then  for  an  isotropic  surface  the  Fourier  transformation  between  the 
roughness  spectral  density  and  the  correlation  function,  as  given  in 
Eq.  (7b)  reduces  to  the  following  Hankel  transform  pair. 


(18a) 

and 


4h2  p(r)  a  2ir  ^  tw(t)JQ(rt)  dt. 


(18b)  w(t)  =  rp(  r)  JG(rt)  dr  . 

wo 

Then  Eq.  (16)  may  be  rewritten 


(19a)  ^  ®S'  =  ~  ^o  k  cos  %  cos  ®s  I  a  yg  I  .1  » 


where 


(19b) 

(19c) 


w(t)  =  2tt 


\  r  p  (r)  J0(rt)dr  , 
<o 


t  =  kQ  sin2  0^-2  sin  8^  sin  0S  cos  4>s  +  sin2 


When  the  average  incoherent  scattering  cross  section  per  unit 
surface  area  is  written  as  in  Eq.  (19a),  the  effect  of  the  surface  rough¬ 
ness  is  contained  entirely  in  the  factor  h2  I.  The  dependence  upon 
incidence  and  scattering  geometries,  polarization  states,  and  surface 
material  properties  is  contained  in  the  factors  cos2  0^  cos2  0S  |ayg[2  • 
As  with  all  objects  whose  size  is  small  compared  to  wavelength,  the 
scattering  cross  section  varies  inversely  as  the  fourth  power  of  wave¬ 
length  . 


Let  us  now  interpret  Eqs .  (19)  physically.  They  indicate  that  the 
average  intensity  of  the  scattered  field  in  a  given  direction  varies  indirect 
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proportion  to  the  surface  roughness  spectral  strength  at  surface  rough¬ 
ness  .frequency,  t  (radians/meter).  The  quantity,  t,  in  turn  is  a  function 
of  the  incidence  and  scattering  directions,  as  given  in  Eq.  (19c).  For 
example,  for  scattering  near  the  specular  direction,  t  ~  0,  the  scatter¬ 
ing  cross  section  is  directly  proportional  to  the  roughness  spectral 
strength  at  and  near  DC.  On  the  other  hand,  the  highest  roughness 
spectral  frequencies  which  can  affect  the  average  incoherent  scatter¬ 
ed  power  are  tnose  at  t  =  2k0.  This  occurs  where  9i  =  0S  =  ir/ 2  and 
<)>s  =  tt,  i.e.,  for  backscattering  near  grazing.  Hence,  the  only  possible 
surface  roughness  spectral  frequencies  which  affect  the  average  scatter¬ 
ed  power  in  any  arbitrary  direction  occur  in  the  range  0  <  t  <  2k0;  any 
higher  roughness  spectral  frequencies  present  in  the  surface  can  have 
no  effect  on  the  scattered  power,  at  least  to  the  first  order  perturbation 
theory  used  here  and  subject  to  the  restrictions  listed  above.  The 
physical  interpretation  of  the  scattering  process  is  similar  to  that  for 
the  sinusoidal  surface  mentioned  previously. 

Let  us  now  examine  two  familar  choices  for  the  surface  rough¬ 
ness  correlation  coefficient,  p(r),  and  find  the  resulting  spectral 
densities  w(t)  and  the  corresponding  quantity  I. 

(1)  Gaussian  Surface  Height  Correlation  Coefficient: 
p(r)  =  exp  (-r2/i2  )  . 

Here,  I  is  referred  to  as  the  surface  height  correlation  length. 

The  quantity  I  for  this  case  becomes 

(20)  I  =  iri2exp(-t2i2/4) 

=  uf2  exp  [ -(kQ2i  2/4)  x  (sin2  0^  -  Zsin  0^  sin  0S  cos  ips  +  sin2  0S)J  . 

Surfaces  having  such  a  Gaussian  surface  height  correlation  coefficient 
are  smoothly  curving  with  derivatives  at  all  points.  The  total  mean 
square  slope  of  the  surface  at  any  point  with  the  above  correlation 
coefficient  is 


(2)  Exponential  Surface  Height  Correlation  Coefficient: 
p(r)  =  exp  (-  |  r|  /i  )  . 
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The  quantity  I  in  this  case  becomes: 


2ttjI  2  2Tri2 

(22)  1= -  = - 

[  1  +  t2!2]5/2  [l  +  k02i2(sin2  8i-2  Sin0iSin0s  +  sin2  0  8)JS/2 

Surfaces  with  exponential  correlation  coefficients  are  jagged  and 
have  many  vertical  facets.  The  mean  surface  slopes  and  all  higher  sur¬ 
face  derivatives  are  infinite,  obviously  because  of  these  vertical  facets. 
As  an  example,  an  urban  area  having  buildings  and  houses  comprising 
its  scattering  surface  has  a  surface  height  correlation  coefficient  which 
behaves  like  the  exponential,  especially  near  the  origin. 

It  should  be  noted  that  for  2k0i  <1,  the  behavior  of  I  for  the  two 
correlation  coefficients  is  not  significantly  different,  except  for  a 
factor  of  2.  Nor  are  the  above  two  correlation  coefficient  models 
necessarily  the  mo,st  applicable  in  all  cases.  Peake  (I960)  found  that 
a  model  which  empirically  fitted  measured  points  for  the  correlation 
coefficient  of  an  asphalt  surface  was 

p(r)  =  [  1  +  20  rzflz 


where  r  is  in  cm.  By  the  Hankel  transformation  between  I  and  p(r), 
in  this  case  I  is  an  exponential  function  of  t. 


C .  Results  for  Vertical  and  Horizontal 
Polarization  States 


We  obtained  above  the  average  incoherent  scattering  cross  section 
for  a  horizontally  polarized  incident  plane  wave  and  for  the  horizontally 
and  vertically  polarized  scattered  components.  We  left  them  in  terms 
of  scattering  matrix  elements,  and  ctvh,  which  are  determined  by 

Eqs .  (9)  and  (10).  The  resulting  cross  sections  and  scattering  matrix 
elements,  ot^y  and  avv,  for  a  vertically  polarized  incident  plane  wave 
may  be  found  immediately  by  noting  that  a  vertically  polarized  plane 
wave  (vertical  or  $  directed  E-field)  is  identical  to  a  horizontally 
directed  H-field  vector.  Hence  the  quantities  othv  and  Qvv  may  be 
obtained  from  avh  and  a  hh  by  interchanging  the  roles  of  the  relative 
constitutive  constants,  «r  and  pr,  and  employing  care  in  observing 
the  resulting  signs  of  the  E-field  vectors  from  a  knowledge  of  the  H- 
field  vectors.  (The  right  subscript  represents  the  polarization  state 
of  the  incident  wave  and  the  left  subscript  represents  the  desired  scat¬ 
tered  wave . ) 
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1  .  Bistatic  scattering  matrix  elements 
for  homogeneous  surface 


(23a)  ahh=- 


(|xr -  I)(fjtrsin%sin0s-cos4>s  >«r|JLr-sin2ei>J«r|jLr-sin20s)  +  ^r2(«r-l)coi 


[  prcosG^  4^«rpr- sin2  [  prcos0s  4 

■  J  «rpr-sin2es] 

^rt^r-^  ^  *rlAr" sin2  -  H-r^r-^)  ^ 

«rpr-sin2  0S 

[prcos  6i  +  j€r|ir-sin26iJ  [  «rcos  0S  +  Jtrpj.- sin2  0S  ] 

M-r(«r-1)  \ 

tj-pr-sin2^  -  «r(pr-l)  > 

«rpr 

-sin2  0s 

[  «rcos6i  + 

«rpr-  sin2  %J  [  prcos  0S  + 

«rpr-sin20s  ] 

(«r-  1  )(trsin0isin0s-cos<}>s  v 

€rpr-sin20i  > 

€rpr-sin2Gs)  +€r2(pr_1)co8 

f.crcos0i  4  rpr- sin2  Gj]  [ 

€rcos0s  4-^€rpr-sin2  0g] 

In  the  above  expressions,  er  and  pr>  the  relative  surface  material 
permittivity  and  permeability,  may  each  be  either  real  or  complex 
(indicating  lossless  or  lossy  material),  and  be  either  greater,  equal  to, 
or  less  than  unity  in  magnitude. 


2 .  Backscattering  matrix  elements  for 
homogeneous  surface 


The  above  elements  are  readily  reduced  to  the  case  of  back- 
scattering  where  4»s  =  tt  and  0S  =  ®i’ 


(24a) 

(24b) 

and 


(Hr' 1 )  t^r-1)  sin2  d{  +  cr|xrJ  -  pr2  («r-l) 

ct  hh  .  , - — — .  2 

[  pr  cos  +  j€rpr  -  sin2  GJ 

Q  vh  =  °  hv  =  ®  > 


(24c) 


(«r_l)  [(*r-^  sin<i  ^i  ^rM’r]  "  €r2  (^r-*) 

a  vv  -  ■  - - 

[«r  cos  Q{  4  J  «rHr  -  sin2  GjJ 
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Bistatic  scattering  matrix  elements 
for  perfectly  conducting  surface 


The  scattering  matrix  elements  for  a  perfectly  conducting  surface 
may  be  readily  determined  from  Eqs .  (23)  by  putting  pr  =  1  and  «r  -*■  «>;* 
this  is  evident  from  consideration  of  the  basic  boundary  conditions  for 


the  two 

cases . 

(25a) 

ahh  =  "  cos  * 

(25b) 

sin  4>s 

a  vh  H  * 

COS  Dg 

(25c) 

sin  4>s 

Q  hv  a  > 

cos 

(25d) 

sin  9^  sin  0g 
°vv  cos  6i  cos 

4  .  Backscattering  matrix  elements  for 
perfectly  conducting  surface 

Here  again,  we  set  4>s  =  ir  and  9S  =  in  the  above  equations. 


(26a) 

ahh  ~  *  > 

(26b) 

avh  =  ahv  =  0  * 

1  +  sin2  0i 

(26c) 

QW-  cos2  9; 

$ 

Actually,  a  perfectly  conducting  surface  implies  that  the  imaginary 
part  of  «r  approaches  infinity;  the  same  result  is  obtained  for  either 
limit,  however. 
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The  average  incoherent  scattering  cross  section  per  unit  surface 
area  is  now  obtained  using  Eqs .  (19)  along  with  the  proper  scattering 
matrix  elements,  a-yg,  for  the  desired  combination  of  vertical  and 
horizontal  polarization  states.  These  are  given  above.  The  factor  I 
used  in  Eq.  (19a)  is  given  in  terms  of  the  surface  roughness  spectral 
density,  as  discussed  previously  and  shown  in  Eq.  (19b). 

Two  outstanding  facts  are  in  evidence  from  the  atove  equations 
for  the  important  case  of  radar  backscattering .  For  one  thing,  the 
cross  polarized  components  are  zero.  It  should  be  pointed  out,  how¬ 
ever,  that  these  components,  as  indicated  above,  are  zero  only  to  the 
first  order.  The  next  order  neglected  terms  in  0-°^  and  are  pro¬ 
portional  to  the  mean  square,  height  and  slope,  kQ2h2  and  s2,  with 
respect  to  and  cr^v.  An  expression  for  these  terms  is  derived  by 
Valenzuela  (1967).  The  lowest  order  terms  retained  here  vary  with 
k04  while  the  next  order  terms  derived  by  Valenzuela  are  of  order  kQ4 
in  frequency.  Thus  our  results  indeed  represent  the  solution  in  the  low 
frequency  limit.  The  second  fact  is  that  backscatter  differs  for  the 
vertical  and  horizontal  states.  The  vertical  component  is  always 
greater  than  the  horizontal  for  |«r|  >  furl*  except  at  normal  incidence 

where  they  are  equal  (i.e.,  Sj  =  0).  This  fact  has  been  observed  experi¬ 
mentally,  but  previous  theories  based  upon  the  less  accurate  tangent 
plane  approximation  have  failed  to  show  any  such  difference. 

(Davies,  ibid.). 

Plots  showing  the  average  incoherent  backscattering  cross  section 
per  unit  surface  area  are  presented  in  Fig.  2  for  perfectly  conducting 
surfaces  along  with  dielectric  surfaces.  The  cross  sections  for  both 
vertical  and  horizontal  polarization  states  are  given.  These  cross 
sections  are  normalized  by  dividing  by  ko2h2  .  The  Gaussian  correlation 
coefficient  model  given  above  was  used  with  two  relative  correlation 
lengths  (i.e.,  koi  =  1.0  and  kQi  =  5.0).  The  results  for  kD|  <1.0  are 
not  significantly  different  in  shape  from  those  for  kQ!  =  1.0,  but  the 
magnitude  of  the  cross  section  itself  varies  in  direct  proportion  to  k02/2, 
as  seen  from  Eqs.  (19a),  (20),  and  (22).  Note  that  as  correlation  length 
(i.e.,  k0f )  increases,  more  of  the  backscattered  power  is  concentrated 
near  the  specular  direction  at  normal  incidence  ar.d  less  near  grazing. 
This  is  physically  reasonable,  because  surfaces  with  a  given  mean 
roughness  height,  h2 ,  are  smoother  when  the  roughness  height  corre¬ 
lation  length,  i,  is  longer.  This  can  be  seen  from  Eq.  (21)  for  the 
mean  square  slope  of  a  rough  surface. 
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'•  vix; 


From  wnufli  ,«fw  ir  .rw  ,*.i.co  *g«« 

(3)  ,’4) 

Fig.  2.  Average  incoherent  backscattering  cross  sections  per  unit 
area  with  linear  polarization  states  from  slightly  rough 
surface  model  vs.  incidence  angle  for  Gaussian  surface 
height  correlation  coefficients,  normalized  to  ko2h2. 

(1)  koi  =  1,  vertical  polarization  states,  (2)  ko i  =  1 .0, 
horizontal  polarization  states,  (3)  koi  =  5. 0,  vertical 
polarization  states,  (4)  koi  =5.0,  horizontal  polarization 
states . 
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D.  Results  for  Circular  Polarization  States 


It  is  relatively  easy  to  transform  from  the  vertical  and  horizontal 
polarization  states  to  the  circular  polarization  states  by  using  standard 
matrix  methods.  As  mentioned  previously,  the  elements  avv,  ahh* 
avh,  and  cr^v  are  directly  proportional  to  scattering  matrix  elements, 
preserving  both  proper  relative  amplitudes  and  phase  differences. 
Hence,  a  set  of  analogous  elements  applicable  for  circular  polarization 
are  determined  and  given  below  for  use  in  Eq.  (19a). 


(27a) 


1 .  Bistatic  scattering  matrix  elements* 


aLR  = 


ahh  +  avv  +  i(ahv  -  avh) 


(27b) 


ahh  "  avv  +  i(a  hv  a  vh) 
aRR - ; - 


(27c) 


„  _  ahh  +  a  vv  “  Uahv  *  avh^ 

aRL - 5 - 


(27d) 


_  ahh  ~  avv  ~  Hahv  +  avh^  • 
ULL  ” 


The  quantities  on  the  right  side  of  these  equations  are  given  in 
Eqs.  (23)  for  a  homogeneous  surface  material  or  by  Eqs .  (25)  for  a 
perfectly  conducting  surface. 


The  right  subscript  represents  the  polarization  state  of  the  incident 
wave  and  the  left  subscript  represents  that  of  the  desired  scattered 
wave . 
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2  .  Backscattering  matrix  elements 


In  the  case  of  backscattering,  it  was  shown  in  Eqs .  (24)  that 
avh  =  ahv  =  0.  Hence  the  above  elements  reduce  to  the  following: 


(28a) 


_  ahh  +  ^vv 
aLR  =  aRL - 1 - 


(28b) 


aLL  =  aRR 


ahh  "  avv 
2 


In  this  case,  the  elements  on  the  right  side  are  given  by  Eqs.  (24)  for 
a  homogeneous  surface  material  and  by  Eq.  (26)  for  a  perfectly  con¬ 
ducting  surface. 

Physically,  cirl  and  aLR  represent  the  "polarized"  components 
in  the  scattered  field;  for  backscattering  from  a  smooth  surface  or 
symmetric  object,  the  scattered  circular  state  is  always  opposite 
sense  from  that  of  the  incident  circular  state.  "Depolarization",  then, 
appears  in  the  elements  oll  and  ctrr,  representing  the  scattered 
circular  state  of  the  same  sense  as  the  incident  state.  It  is  important 
to  note  that  depolarization  occurs  for  the  circular  states  for  back- 
scattering  even  though  it  does  not  occur  for  the  linear  states  (i.e., 
ahv  =  avh  =  0):  this  is  due  to  the  non- zero  difference  between 
and  flyv  • 

Curves  showing  the  normalized  average  backscattering  cross 
section  per  unit  surface  area  as  a  function  of  incidence  angle  are 
presented  in  Fig.  3.  These  represent  perfectly  conducting  and  dielectric 
surfaces  for  both  the  "polarized"  and  "depolarized"  components .  The 
Gaussian  surface  height  correlation  coefficient  model  is  used  with  two 
values  of  the  correlation  length,  kof  =  1.0  and  5.0.  It  must  be 
emphasized  again  that  all  of  these  curves  were  derived  from  a  first- 
order  perturbation  theory;  higher  order  terms  in  k0h  and  s  (relative 
rms  roughness  height  and  slope)  have  been  neglected. 
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(3)  (4) 

Fig.  3.  Average  incoherent  backscattering  cross  sections  per 
unit  area  with  circular  polarization  states  from  slightly 
rough  surface  model  vs.  incidence  angle  for  Gaussian 
surface  height  correlation  coefficient,  normalized  to 
kQ  h  .  (1)  koi  =  1.0,  opposite  sense  polarization  states, 

(2)  koi  =  1.0,  same  sense  polarization  states, 

(3)  koi  =  5.0,  opposite  sense  polarization  states, 

(4)  kQi  =  5.0,  same  sense  polarization  states. 
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E •  Results  for  Arbitrary  Linear 
Polarization  States 


In  many  experimental  situations,  the  transmitting  and  receiving 
antennas  employ  linear  polarization  states,  but  they  may  not  necessarily 
be  oriented  along  either  the  vertical  or  the  horizontal  polarization 
directions  defined  by  the  surface  coordinate  system.  Define  as  the 
angle  between  the  unit  vector  and  the  desired  incident  E-field  com¬ 
ponent  in  the  direction  toward  the  unit  vector  Define  t}s  in  the  same 

way  with  respect  to  0S  and  <j>s  Then  the  matrix  element  a  for 
bistatic  scattering  can  be  written  s  * 


(29)  a^.  =  avv  cos  cos  +  Qvh  sin  ’H  cos 

+  ahv  cos  ’’H  s^n  ^s  +  a^h  sin  sin  tis 


1 .  Backscattering  with  aligned  and 
crossed  linear  polarized  antenna 

Often  the  transmitter  and  receiver  antennas  are  aligned  so  that 
Tis  =  TH -  For  the  receiver  antenna  polarization  direction  orthogonal  to 
this,  qg  =  t]{  +  it/Z.  Employing  these  angles  and  also  the  fact  that 
avh  =  ahv  =  0  for  backscattering,  we  arrive  at  the  following  two  elements. 

(30a)  “mni  =  QVV  cos2  "Hi  +  ahh  sin2  , 

and 

(30b)  am+ir/2,ili  ='(avv  -  ahh>  sin  ’li  co8  ^  * 


These  elements  can  be  used  in  Eq.  (19a)  to  determine  the  back- 
scattering  cross  section.  The  elements  ahh  and  O^v  to  be  used  in  the 
right  side  of  the  above  are  given  in  Eqs  .  (24)  for  a  homogeneous 
surface  material  and  Eq.  (26)  for  a  perfectly  conducting  surface. 
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2  .  Angularly  averaged  backscattering 
with  aligned  and  crossed  linear 
polarized  antennas 


In  many  cases  aligned  and  crossed  linear  transmitting  and 
receiving  antennas  may  be  employed,  but  the  antenna  angle,  t^,  with 
respect  to  the  surface  vertical  may  be  rotating  and  unknown  at  a  given 
instant  of  time.  This  may  arise  in  the  case  of  satellite  observations 
of  a  planetary  surface.  In  other  situations  it  can  result  from  uncon¬ 
trollable  Faraday  rotation  of  polarization  due  to  the  ionosphere.  An¬ 
other  situation  where  an  average  quantity  is  desired  is  the  case  where 
a  relatively  short  pulse  of  linear  polarization  sweeps  past  a  planetary 
surface,  illuminating  an  annular  portion  of  the  sphere.  The  angle, 
between  the  linear  direction  and  the  local  plane  of  incidence  at  a  given 
point  on  the  annulus  varies  continuously  between  0  and  2tt  in  this  case. 

In  the  above  cases,  the  quantities  desired  are  [ ot-rj- f 2  and 
laTli+ir/2,  I  2  ’  averaged  by  integrating  over  r|i  from  0  to  2-rr  and 
dividing  by  2tt  .  These  quantities  become 


(31a)  <|aTliT1.|2>  =  -i[3|avv|2  +  2Re(avvahh)  +  3|ahh|2  J  . 

and 


(31b) 


■  I  a  I  avv  ■  ahh  I 

<!aT1i+7r/ 2, nil  >=  - g - 


The  notation  Re(x)  denotes  the  real  part  of  x. 


F .  Comparison  with  Measured  Results 


It  is  difficult  to  find  measured  data  with  which  to  compare  the 
preceding  theoretical  results  for  a  slightly  rough  surface  for  two  reasons, 
(i)  It  is  necessary  that  roughness  height  be  everywhere  small  compared 
to  wavelength,  (ii)  It  is  difficult  experimentally  to  separate  the  average 
coherent  scattered  power  from  the  incoherent;  in  most  cases,  no  attempt 
is  made  to  do  this.  For  backscattering,  the  coherent  contribution  to 
the  total  power  dominates  near  normal  incidence.  The  results  p~esented 
in  this  section  apply  only  to  the  incoherent  scattered  power. 

Recently,  average  backscattering  cross  section  measurements 
were  reported  by  Wright  (1966)  from  a  fresh  water  surface  at  X-band. 
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Capillary  waves  v/ere  generated  under  controlled  conditions  such  that 
the  wave  amplitudes  were  small  compared  to  radar  wavelength.  Return 
for  the  both  vertical  and  horizontal  states  were  obtained,  and  are  shown 
in  Fig.  4  as  a  function  of  incidence  angle.  Along  with  these  measured 
curves,  theoretical  curves  are  drawn,  using  the  equations  derived  above 
and  the  parameters  of  the  water  surface  estimated  by  Wright.  These  are 
«r  -  55(1  +  i.55),  kQh  =  .05,  kQ i  =  2.  It  should  be  noted  that  the  measured 
results  were  presented  not  in  terms  of  backscattering  cross  section  per 
unit  area,  but  backscattering  cross  section  in  square  meters.  Since  the 
area  of  the  illuminated  surface  was  not  given,  the  absolute  value  of  cr°v 
and  irgh  cannot  be  ascertained.  Hence,  the  magnitudes  shown  in  the 
figure  are  made  to  conform  to  the  theoretical.  The  significant  observa¬ 
tion  is  the  agreement  in  functional  dependence  between  the  measured 
and  theoretical  curves. 

Another  set  of  measurements  on  an  asphalt  surface  at  X-  and 
Ka-bands  allow  another  comparison  with  the  model.  (See  Fig.  5). 

Actual  measurements  of  surface  statisitics  and  dielectric  constant 
were  made  in  this  case  and  used  to  compute  the  curves  for  the  slightly 
rough  surface  model,  (Cosgriff,  I960).  Again  the  functional  dependence 
on  frequency,  angle,  and  polarization  is  in  good  agreement.  Even  the 
somewhat  poorer  agreement  in  absolute  value  is  reasonable  in  view  of 
the  large  number  of  surface  parameters  which  had  to  be  estimated 
from  rather  small  samples. 


III.  VERY  ROUGH  PLANAR  SURFACE 
A.  Introduction 


The  remaining  class  of  rough  surfaces  which  can  be  treated 
analytically  is  the  one  whose  memoers  have  a  large  scale  roughness 
such  that  the  local  surface  radii  of  curvature  over  nearly  all  surface 
regions  are  much  greater  than  wavelength.  In  most  cases  of  practical 
interest,  the  rms  roughness  height,  h,  will  then  also  be  larger  than  a 
wavelength.  Then  all  of  the  scattered  power  from  such  a  surface  is 
incoherent,  (i.e.,  the  phase  angle  of  the  scattered  field  becomes 
uniformly  distributed  between  0  and  2tt,  and  hence  the  average  value  of 
the  scattered  field  is  zero;  only  the  average  scattered  power  is  non¬ 
zero.) 

When  this  radius  of  curvature  criterion  is  satisfied,  the  scatter¬ 
ing  problem  can  be  solved  by  several  optics  techniques.  If  these  are 
used  properly,  they  all  lead  to  the  same  result.  The  value  of  examining 
more  than  one  approach  to  the  problem  lies  in  the  physical  insight  into 
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Fig.  4. 


Measured  average  back- 
scattering  cross  section 
per  unit  area  from  slightly 
rough  fresh  water  surface 
at  X-band  vs.  incidence 
angle  for  vertical  and  hori¬ 
zontal  polarization  states. 
Dashed  curves  represent 
slightly  rough  surface  model 
using  parameters  measured 
from  fresh  water  surface. 


ANGLE  FROM  VERTICAL  (  DEGREES)  ANGLE  FROM  VERTICAL  I  DEGREES ) 


Fig.  5.  Measured  average  backscatte ring  cross  section  per 
unit  area  from  slightly  rough  asphalt  surface  at  X- 
and  Ka-bands  vs.  incidence  angle  for  linear  polari¬ 
zation  states.  Dashed  curve  is  slightly  rough  sur¬ 
face  model  using  measured  parameters  p(r)  =  (l+20r2)-3/2 
where  r  is  in  cm;  kQh  -  0 . 085  at  X-band  and  0.30  at 
Ka-band;k04  =2.  0-48  at  X-band  and  1.7  at  Kg-band; 
«r~4.3  +  i0.1  at  X-band  and  2.5  +  iO  .65  at  Ka  -  band. 
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the  scattering  process  which  each  provides.  At  present,  three  different 
optics  approaches  have  appeared  in  the  literature.  Rather  than  retrace 
a  derivation  in  detail  here,  a  critique  of  these  existing  analyses  will  be 
made,  along  with  the  physical  explanation  of  each.  Finally,  the  common 
results  of  these  analyses  will  be  presented  and  curves  for  two  models 
will  be  shown  for  comparison. 

All  of  the  three  optics  approaches  apply  to  the  class  of  rough 
surfaces  which  satisfy  (at  least)  the  following  restrictions: 

(a)  The  local  surface  radii  of  curvature  at  nearly  every  point 
on  the  surface  are  considerably  larger  than  wavelength. 

(b)  3£/8x,  8f,/8y  <1.0,  i.e.,  surface  slopes  are  relatively 
small.  Only  when  this  is  t.-ue  can  one  neglect  shadowing 
and  multiple  scattering  over  most  of  the  range  of  incidence 
and  scattering  angles  . 

(c)  <(8^/8x)2>  =  <(9^/8y)2>,  i.e.,  the  roughness  is  isotropic. 

The  restriction  is  not  necessary  to  obtain  an  answer,  but  is 
employed  in  this  paper  for  simplicity. 

(d)  L  >  >  i  >  >  X,  i.e.,  the  dimensions  of  the  illuminated  area 
are  much  greater  than  the  roughness  correlation  length,  i, 
which  in  turn  is  much  greater  than  wavelength.  The  second 
portion  of  the  inequality  is  automatically  fulfilled  if  (a)  is 
true.  The  first  part  of  the  inequality  precludes  the  possibility 
of  looking  at  only  a  part  of  a  hill  or  wave,  but  rather  it  de¬ 
mands  that  several  hills  or  waves  be  included  in  the  illuminated 
area. 

(e)  k02h2  >>  1,  i.e.,  mean  square  roughness  height,  h2,  is  of 
the  order  of,  or  greater  than,  wavelength  squared. 

(f)  The  surface  height  correlation  coefficient,  p(r),  must  be 
parabolic  in  behavior  at  the  origin,  i.e.,  p(r)-  j)l  -  r 2  /  i2  , 
where  r  is  the  separation  between  two  points  on  the  surface 
and!  is  termed  here  the  correlation  length. 

The  first  two  restrictions  above  permit  use  of  the  tangent  plane 
(or  physical  optics)  approximation.  This  approximation  involves  writing 
the  scattered  field  at  the  surface  as  the  product  of  the  incident  field  and 
the  proper  reflection  coefficient.  Restriction  (f)  is  imposed  in  order  to 
avoid  mathematical  difficulties  which  lead  to  erroneous  results  when 
one  employs  a  correlation  coefficient  whose  behavior  is  linear  at  the 
origin,  (Barrick,  1966).  The  second,  third  and  fourth  restrictions  are 
identical  to  those  given  previously  for  the  slightly  rough  surface. 


28 


Although  the  slightly  rough  surface  analysis  of  the  preceding  part 
takes  into  account  multiple  scattering  and  shadowing,  the  normal  optics 
techniques  used  for  very  rough  surfaces  ignore  these  multiple  inter¬ 
action  effects.  Generally,  they  are  higher-order  effects  so  long  as  the 
surface  slopes  are  not  too  precipitous  and  the  scattering  angles  are  not 
too  close  to  grazing.  No  known  correction  factor  for  multiple  scattering 
has  been  computed,  but  shadowing  corrections  have  been  formulated 
recently.  Beckmann  (1965)  first  formulated  the  "shadowing  function." 
This  was  later  revised  by  Brockelman  and  Hagfors  (1966)  and  Wagner 
(1967).  The  latter  shows  that  the  correction  factor  is  less  than  3  dB 
for  angles  of  incidence  and  scattering  more  than  10°  above  grazing  and 
an  rms  total  surface  slope,  tan_1s  =  tan_1(2h/i),  of  less  than  27°.  In 
view  of  the  restrictions  placed  upon  the  class  of  surfaces  considered 
here,  shadowing  corrections  will  not  be  introduced  in  this  paper. 


B .  Review  and  Interpretation 
1 .  Physical  optics  approach 

The  first  approach  to  appear  historically  is  also  probably  the  most 
mathematically  thorough  and  satisfying.  By  the  same  token,  it  probably 
sheds  the  least  physical  insight  on  the  scatter  mechanism  because  of  the 
mathematical  detail.  This  method  employs  a  physical  optics  integral 
expression  for  the  scattered  field.  The  total  fields  at  the  surface 
appearing  in  the  integrand  are  determined  from  the  tangent  plane  approxi¬ 
mation  mentioned  above.  The  first  comprehensive  treatment  of  this 
problem  for  a  perfectly  conducting  surface  was  made  by  Isakovich  ( 1952) 
who  employed  a  vector  physical  optics  formulation.  This  work  appears 
to  have  been  overlooked  by  Western  investigators,  for  Davies  (1954) 
solves  the  same  problem  by  the  same  techniques,  but  employs  a  less 
complete  scalar  formulation.  Beginning  about  I960,  a  variety  of 
Western  and  Soviet  investigators  began  to  extend  these  results.  Among 
these  are  Hagfors  (I960,  1964),  Beckmann  (1963),  Hughes  (1962),  Fung 
(1964),  Hay  re  (1961),  Daniels  (196  1),  Semenov  (1965),  and  Stogryn  (1967) 
Semenov  appears  to  have  been  the  first  to  have  solved  the  problem  for 
scattering  from  a  surface  of  homogeneous  material,  cr,|ir,  using  a  vector 
formulation,  which  accounts  properly  for  polarization.  The  results  of 
Semenov  have  been  derived  independently  and  simultaneously  by  Barrick, 
and  the  two  are  in  agreement,  after  some  algebraic  simplification  of 
Semenov's  results.  Stogryn's  solution  is  the  same  as  that  of  Semenov 
(although  his  notation  is  more  complicated),  and  hence  serves  as  a  third 
check  on  Semenov's  analysis. 

Semenov  expresses  the  scattered  vector  field  (for  both  vertical 
and  horizontal  incident  polarization  directions)  in  terms  of  a  physical 
optics  integral.  For  example,  for  an  incident  plane  wave  whose  E-field 
is  horizontally  polarized,  the  integral,  after  application  of  the  tangent 
plane  approximation,  can  be  written  as  follows:  (See  Fig.  1). 
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where 


RQ  =  distance  from  origin  to  observation  point, 

r  =  xx  +  y  y  +  £(x, y)z  =  distance  from  origin  to  point 

on  rough  surface, 

k{,  icg  -  unit  constant  vectors  pointing  in  direction  of 
incidence  and  scattering, 


Cx»  Cy  “  local  surface  slopes  in  x  and  y  directions  at 
surface  point  t,(x,  y),  i.e.,  8£/8x  and  8  £/8y  • 

The  factor  F(£x>  £y)  is  a  function  of  the  local  normal  to  the  surface 
and  of  the  local  Fresnel  reflection  coefficients  at  each  surface  point.  It 
can  be  expressed  as  follows: 


(33a)  F(Cx.£y)  = 


jl  +  !»XZ  +  £y2 

sin2  Vi 


.  [y.$i  X  n)  Gx  +  y.  (ki  cos  y{  +  n)GJ  , 


where 


(33b)  Gx  =  [cos  Yi  -  cos  y8  RjYi)]  [  (^  X  n)  -  £s  .(iq  X  n)£s] 

A  A  A 

+  [  1  +  Rx(Yt)J  l^s  X  (^  +  n  cos  Yi). 


(33c)  G„  =  £1  +R„(Yi>]  £$> 
+  [cos  Yi  -  cos  Ys 


+  n  cos  Yi)  -  £s*^i  +  “  cos  Yi)^sJ 
R„  <Yi>]  t^i  cos  y8  -  (^s'fy  * 


(33d)  n  = 


A  _  A  A 

-^c*  -  4vy  +  z 


=  local  unit  normal  vector  to  surface  at 
point  £(x,  y), 


Yi*  Ys  "  l°cal  incidence  and  scattering  angles  in  planes  of  incidence 
and  scattering  with  respect  to  local  surface  normal,  (i.e., 
cos  Yi  =  ■  ^i*n,  cos  Ys  =  ^s*”)» 
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Rx(Vi)«  R„  (Yi)  “  Fresnel  reflection  coefficients  at  local  surface 
point  £(x,y)  with  incidence  angle,  Yi ' 

As  one  can  observe,  the  integrand  consists  of  the  factor  F(£x.£y), 
which  varies  with  the  local  surface  slopes  from  point  to  point  on  the 
surface,  and  an  exponential  phase  factor  containing  the  phase  difference 
between  the  local  surface  point  and  the  scattered  field  point.  The  first 
factor  is  to  be  removed  from  the  integrand  as  a  constant;  the  justification 
for  this  is  the  stationary  phase  (or  specular  point)  argument.  This 
argument  maintains  that  for  large  kQ  in  the  exponential  factor,  the  only 
surface  regions  which  contribute  to  the  integral  are  those  smoothly 
curving  portions  which  are  in  a  position  to  specularly  reflect  into  the 
^esired  scattering  direction.  The  direction  of  the  surface  normal, 
nSp,  and  the  local  incidence  angle,  Yi  =  Ys  =  1  •  at  these  specular  points 
are  readily  determined  by  bisecting  the  angle  between  the  incidence 
propagation  direction  and  the  desired  scattering  propagation  direction. 

The  constant  factor,  F(t,x  ,t,yg  ),  is  then  evaluated  using  the  slopes 
of  the  surface  at  these  specular  {joints  . 

The  remaining  integral  containing  the  exponential  phase  factor 
is  not  solved  by  the  same  specular  point  or  stationary  phase  approach, 
however;  the  scattered  field  is  squared  to  give  scattered  power  or 
intensity,  and  then  it  is  averaged.  Averaging  consists  of  multiplying 
the  double  integral  by  the  joint  probability  density  functions  for  the 
surface  height  at  two  different  surface  points  and  then  integrating  over 
these  two  random  variables.  The  order  of  integration  is  interchanged, 
and  the  integrand  of  the  former  double  integral  is  averaged  first.  This 
interchange  is  justified  mathematically  (See  Middleton,  I960)  so  long  as 
restriction  (f)  is  applied.  With  the  aid  of  restriction  (e),  the  remaining 
integral  is  solved  in  a  straightforward  manner;  nearly  all  of  the  above 
mentioned  references  carry  out  this  latter  process  correctly. 

The  step  of  averaging  under  the  integral  sign,  though  mathe¬ 
matically  correct,  obscures  the  physical  understanding  of  the  scattering 
mechanism.  It  is  also  somewhat  less  than  satisfying  to  justify  the 
removal  of  a  complicated  factor  in  the  integrand  by  the  specular  point 
or  stationary  phase  argument,  but  then  never  to  actually  employ  this 
stationary  phase  principle  to  solve  the  remaining  integral.  However, 
the  two  approaches  to  be  discussed  subsequently  do  justify  the  validity 
of  this  procedure  by  actually  showing  that  the  specular  point  theories 
give  exactly  the  same  result,  at  least  for  a  Gaussian  distributed  surface 
height. 

No  one  has  yet  justified  correctly  a  more  exact  method  of  evaluating 
this  integral  than  the  stationary  phase  approach.  This  is  obviously  due  to 
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the  extremely  complicated  form  of  the  factor  F(£,x»Cy)*  At  least  two 
previous  analyses  have  proposed  alternate  schemes  for  evaluating  the 
integral.  Barrick  (1965)  suggested  averaging  this  factor  independent 
of  the  exponential  factor;  hence,  F(£x,  £y)_would  be  removed  as  a  con¬ 
stant  evaluated  at  the  mean  slopes,  i.e.,  F(0,  0).  Fung  (1965)  proposed 
expanding  the  factor  in  a  series  in  £x,  £  y,  and  averaging  each  term 
separately.  One  can  show  that  the  resultant  series  cannot  be  truncated, 
since  higher  order  terms  are  of  the  same  order  as  the  first  terms. 

This  is  true  even  for  vanishingly  small  values  of  surface  slopes,  £x 
and  £y.  Both  of  these  approaches  are  inconsistent  with  the  stationary 
phase  principle,  and  therefore  both  can  have  no  theoretical  justification. 

Recently,  Kivelson  and  Moszkowski  (1965)  have  shown  that  for 
backscattering  one  can  reduce  the  physical  optics  integral  directly  to 
the  probability  density  function  for  the  surface  normal;  this  thus 
provides  a  valuable  connection  to  the  ray  optics  approach  of  the  next 
section.  Gaussian  surface  statistics  are  not  necessary  for  their 
analysis,  but  the  generalized  form  selected  for  their  conditional  sur¬ 
face  height  probability  density  function  is  still  somewhat  restrictive. 

The  joint  probability  density  function  which  they  form  is  then  not 
symmetric  in  the  two  height  random  variables,  £  and  as  it  must 
be,  except  for  the  special  case  of  Gaussian  statistics. 


2 .  Ray  optics  approach 

Recently,  Muhleman  (1964)  brought  an  older  -  but  better  under- 
stood-theory  to  bear  on  the  very  rough  surface  problem.  The  "ray 
optics"  approach  (following  his  terminology)  provides  a  more  readily 
and  simply  explainable  interpretation  of  the  scattering  process,  even 
though  mathematically  it  is  less  exact  and  elegant.  The  rough  surface 
is  initially  approximated  by  a  grid  of  small  flat  planar  elements,  all  of 
which  are  connected  to  form  the  rough  surface.  Each  element  reflects 
incident  power  specularly,  and  the  direction  of  reflection  from  a  given 
facet  is  determined  by  the  direction  of  its  normal.  Hence  the  amount 
of  power  scattered  into  a  given  direction  is  directly  proportional  to  the 
number  of  facets  whose  normals  are  pointed  in  the  proper  direction 
(i.e.,  the  direction  bisecting  the  angle  between  the  incidence  and  scat¬ 
tering  directions)  divided  by  the  total  number  of  facets.  This  quantity 
is  the  probability  density  function  for  the  surface  normal.  It  is 
expressible  in  terms  of  the  probability  density  function  for  the  surface 
slopes,  siijice  the  surface  slopes  are  perpendicular  to  the  local  normal. 

Muhleman  dealt  with  the  problem  of  bistatic  scattering  from  a 
perfectly  reflecting  surface,  but  the  problem  of  the  homogeneous  surface 
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material  is  just  as  easily  solved.  Recall  that  in  the  latter  case,  the 
components  of  the  local  E-field  incident  upon  a  planar  facet  in  and 
perpendicular  to  the  plane  of  incidence  each  have  different  Fresnel 
reflection  coefficients,  R„  and  Rx.  Referring  to  Fig.  1,  we  can 
determine  the  vertical  and  horizontal  reflection  coefficients  for  a 
properly  oriented  *"acet  which  reflects  power  incident  from  direction 
into  direction  Bg.^s’  Vertical  and  horizontal  here  refer  to  the 
polarization  directions  of  incident  and  scattered  fields  with  respect 
to  the  mean  surface,  i.e.,  the  x-y  surface,  and  not  with  respect  to 
the  local  surface  facet.  Let  us  define  Rvfc  as  the  reflection  coefficient 
between  the  horizontally  polarized  incident  (right  subscript)  and  the 
vertically  polarized  reflected  (left  subscript)  E-field  components  with 
mean  incidence  and  scattering  angles  B^,  0S,4>S.  The  other  three  co¬ 
efficients  are  defined  in  the  same  manner.  Then  assuming  that  scat¬ 
tering  in  these  directions  is  produced  by  planar  facets  properly 
oriented  to  specularly  reflect,  these  vertical  and  horizontal  reflection 
coefficients  are  determined  from  straightforward  geometrical  con¬ 
siderations.  They  are: 


(34a) 


Rhh 


sin  8[  sin  B8  sin2  <t>s  Rl(  (i  )  +  az  ag  Rx(i  ) 
4  sin2 i  cos2  i 


(34b) 


Rvh  =  -  sin  <>8 


a 2  sin  8S  R,j(i)  -  a)  sin  0^  Rx(i) 
4  sin2  t  cos2  l 


(34c) 


RhV  =  -  8in  *8 


a2  sin  08  Rx(i)  -  aj  sin  0*  R,,  (t) 
4  sin2 1  cos2 i 


(34d) 


Lvv 


sin  sin  08  sin2  <t>s  Rx  ( ^)  -  aj  sin  0t  R„  (0 
4  sin2  t  cos2  i 


where 


(35a) 

and 

(35b) 


Ru  (Tl)  = 


Ra(ti)  = 


<  r  cos  ^ 

«rHr  -  sin2  Tl 

<r  COS  T|  + 

J«rHr  -  sin2  r) 

Jir  COS  T)  -  , 

«rHr  -  sin2  r| 

Jlr  COS  Tl  +  ^ 

«r^r  ■  sin2  r| 
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and  where  i  is  the  local  angle  of  incidence  at  the  specular  reflecting 
facets,  and  is  given  by 

(36)  cos  i  =  —  {  1  -  sin  8t  sin  0S  cos  t)>s  +  cos  0^  cos  0s)  • 


The  quantities  a.2  and  aj  involve  the  scattering  angles  and  are  given 
subsequently  in  Eqs .  (54).  As  a  check,  the  above  reflection  coefficients 
for  backscattering  reduce  to  the  following: 


(37)  ^•hh’  ^vv  ~  R||(0),  R-x(^)*  ^hv*  ^-vh  ~ ®  • 


Since  polarization  differences  were  not  considered  by  Muhleman, 
the  above  quantities  squared  must  be  inserted  as  factors  in  his  results 
to  obtain  power  scattering  for  the  vertical  and  horizontal  polarization 
directions  as  defined  in  Fig.  1. 

Hagfors  (1966)  showed  that  the  probability  density  function  for 
surface  slopes  appearing  in  Muhleman's  result  for  scattered  power 
can  be  easily  related  to  the  surface  height  probability  density  function 
when  the  surfaces  are  Gaussian.  In  doing  so,  he  showed  that  for 
backscattering,  one  obtains  the  same  result  as  from  the  physical 
optics  analysis  discussed  above.  It  is  readily  shown  that  the  two 
methods  are  equivalent  for  bistatic  scattering  also.  [Note,  it  was 
pointed  out  by  Hagfors  that  due  to  an  oversight,  oae  must  divide 
Muhleman's  result  by  the  cosine  of  the  angie  between  the  local  surface 
normal  and  the  mean  plane  normal,  i.e.,  by 


(38) 


.  A  A  _ 

cos  (n,  z)  = 


cos  8j  4  cos  8S 
2  cos  i 


•  J 


The  ray  optics  approach  sheds  much  needed  insight  into  the 
scattering  process  from  a  rough  surface.  As  a  rigorous  derivation, 
however,  the  method  can  be  seriously  questioned.  Such  an  approach 
considers  only  scattered  energy  or  power,  and  not  the  fields  themselves; 
hence  phase  relationships  between  fields  from  different  regions  are 
ignored  from  the  outset.  The  ray  optics  method  assumes  specular 
reflection  from  a  small  planar  facet  in  the  form  of  a  non-diverging 
tube  or  parallel  rays.  Later  when  one  permits  the  planar  facet  to 
shrink  to  zero,  the  paradox  between  ray  optics  and  modern  geometrical 
optics  theory  becomes  obvious.  Diffraction  effects  and  patterns,  due  to 
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phase  interference,  and  divergence  of  rays  due  to  curved  phase  fronts 
are  ignored.  The  method,  nonetheless,  serves  our  purpose  here  be¬ 
cause,  under  the  restrictions  mentioned  previously,  we  are  dealing 
with  incoherent  scattered  power.  This  means  that  it  is  sufficient  to 
add  power  scattered  from  different  portions  of  the  surface,  rather 
than  the  fields  themselves.  Hence  ray  optics  is  justified  for  this 
problem,  and  provides  needed  insight  into  the  scattering  process. 


3  .  Geometrical  optics  or  stationary 
phase  approach 

The  third  approach  is  the  application  of  the  stationary  phase 
principle  to  the  Kirchoff  integral  for  the  complex  scattered  field. 

The  result  shows  that  scattering  from  a  portion  of  a  quadric  curving 
surface  does  indeed  radiate  specularly.  The  scattering  cross  sections 
of  such  a  curved  surface  is 


(39)  cr  =  tt  |Rj  R2|  , 


where  Rx  R2  is  the  product  of  the  two  principal  radii  of  curvature  at 
the  specular  or  stationary  phase  point.  A  very  rough  surface  con¬ 
sists  of  many  such  specular  points.  Everyone  has  observed  the 
dancing  specular  points  on  a  rippling  lake  surface  on  a  moonlight  night. 
However,  neither  of  the  previous  theories  directly  relate  the  total 
scattered  power  to  the  number  of  specular  points  and  their  radii  of 
curvature,  even  though  the  physical  optics  method  justifies  the 
simplification  of  the  integrand  on  the  basis  that  this  stationary  phase 
principle  is  valid.  It  is  therefore  instructive  to  examine  the  scattering 
problem  formulated  on  this  basis  and  note  that  one  arrives  at  the  same 
result  as  the  preceding  two  approaches.  In  addition,  this  third 
technique  provides  further  insight  by  explicitly  deriving  the  average 
number  of  specular  points  on  a  rough  surface  and  their  average 
Gaussian  curvature,  Ri  Rz  • 

Kodis  (1966)  formulated  the  rough  surface  problem  rigorously 
in  this  manner  and  showed  by  stationary  phase  that  the  scattering  cross 
section  for  a  rough  surface  can  be  expressed,  as  one  intuitively  expects, 
as  follows: 

N 

(40)  <r=  ^  Z  C I  RliRai I  •  I  Ri jR*j  I  >  * 

i,  j=1 
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where  the  total  number  of  specular  points  is  N  and  the  Gaussian  cur¬ 
vature  at  the  i-th  points  is  RiiR2i*  The  quantity  -  4>j  contains  the 
path  length  difference  between  the  i-th  and  j-th  specular  points. 
Employing  the  restriction  p,  h  >  X  ,  one  can  show  that  for  a  very  rough 
surface,  the  specular  points  are  randomly  and  widely  separated  in 
terms  of  wavelength,  and  therefore,  <j>i  and  4>j  are  uniformly  distri¬ 
buted.  Hence  the  average  scattering  cross  section  becomes 

N 

(41)  <tr>  =  tt  <  )  I  RiiRzi  1  >  • 

Lj 

i=l 

Assuming  now  that  one  can  assign  an  average  Gaussian  curvature  to 
the  specular  points,  one  has 


(42)  <<r>  =  tr  <N>  <|Rj  R2  |>  , 


which  one  expects  for  incoherent  scattering.  The  average  scattering 
cross  section  per  unit  area,  <r°g  for  incident  and  scattered  polarization 
states  &  and  y,  can  now  be  expressed 


(43)  <t°b  =  irn  <|RiR2|>  |Ryb|2 

where 


(44) 


is  the  average  number  of  specular  points  per  unit  area  and  Ryg  is  the 
reflection  coefficient  fcr  the  surface  at  a  specular  point.  For  the 
vertical  and  horizontal  polarization  states  (i.e.,  y,  &  =  v,  h),  these 
coefficients  are  given  in  Eqs.  (32). 

BarriCk  (1968)  used  the  above  formulation  and  obtained  expressions 
for  n  and  <|kiR2|>.  These  results  are  identical  to  those  of  the  previous 
two  approaches.  While  confirming  the  general  validity  of  and  similarity 
among  the  various  high  frequency  optics  techniques  applied  rough  sur¬ 
faces  this  third  approach  also  provides  valuable  insight  into  the  scat¬ 
tering  process,  which  is  complementary  to  that  of  the  other  two.  In 
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particular,  one  finds  that  the  average  number  of  specular  points  is 
proportional  to  the  probability  density  function  for  the  surface  slopes 
evaluated  at  the  required  slopes  at  the  specular  point;  this  quantity 
is  greatest  where  the  local  normal  at  the  specular  points  (or  bisector 
between  incidence  and  scattering  directions)  is  near  the  mean  normal, 
or  vertical.  This  is  reasonable  from  consideration  of  the  moonlit  lake 
example.  The  mean  value  of  the  radii  of  curvature  product  at  the 
specular  points,  on  the  other  hand,  becomes  greater  as  the  local 
normal  is  farther  removed  from  the  mean  normal. 


C.  Results 


1 .  General  form 

Since  the  solutions  for  all  three  approaches  are  essentially  the 
same  for  the  Gaussian  surface,  we  shall  present  one  result  here. 

It  is  expressed  in  a  form  which  proceeds  naturally  from  the  physical 
optics  approach  and  is  obtained  from  the  derivation  of  Semenov  (1965). 
The  average  scattering  cross  section  per  unit  surface  area,  for 
incident  polarization  state  6  and  scattered  state  y,  is  written  as  follows: 

<45>  -V  !pvbI2j  • 


The  quantity  Py6  ^as  exactly  the  same  significance  as  oiyg  for  the 
slightly  rough  surface  model:  it  is  directly  proportional  to  the  scatter¬ 
ing  matrix  element  relating  the  incident  field  of  polarization  state  5 
to  the  scattered  field  component  of  polarization  state  y  •  These  quantities 
will  be  given  below  for  the  vertical  and  horizontal  as  well  as  circular 
states . 

The  factor  J  is  related  to  the  roughness  statistics  in  the  following 
manner: 


(46) 


r  JD(tr)  M^;(iu,  -iu;r)  dr. 


The  quantity  t  is  a  function  of  the  scattering  angles  and  is  given  in 
Eq.  (19c).  The  symbol  u,  also  a  function  of  these  angles,  is  defined  as 

(47)  u  =  kD(cos  ©i  +  cos  0B)  . 
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The  function  M^i(iu,  iv;r)  is  the  characteristic  function  of  the  surface 
height  random  variables  £  and  measured  at  surface  points  (x,  y) 
and(x',y').  (Formally,  it  is  the  double  Fourier  transform  of  the  joint 
probability  density  function,  P(£,V;r);  r  is  the  horizontal  distance 
between  the  surface  points; 

(48)  r  =,  (x-x1)2  +  (y-y1)2  . 


Therefore,  in  order  to  obtain  a  closed  form  result  for  the  scattering 
cross  section,  one  must  either  know  or  assume  a  form  for  the  surface 
height  joint  probability  density  function.  (For  the  slightly  rough  sur¬ 
face,  recall,  one  has  to  choose  a  form  for  the  surface  height  correlation 
coefficient.)  We  shall  select  two  common  probability  density  function 
models  for  the  sake  of  comparison  and  find  the  quantity  J. 


a.  Gaussian  surface  height 
joint  probability  density 

Using  the  standard  Gaussian  joint  probability  density  function 


(49) 


PG(t.t')  = 


2ir  h2[  1  -  p2(r)]2 


[C  2  £'p(r)+  £|2J 

2h2[  1-p2  (r)J 


along  with  restriction  (f)  on  the  form  of  the  form  of  the  correlation 
coefficient  p(r)  near  the  origin  and  restriction  (e)  on  the  size  of  the 
mean  square  roughness  height,  h2 ,  compared  to  wavelength,  we  can 
readily  solve  Eq.  (46)  for  J. 

J_  t2 

(50)  J  =  — — —  e’  **  , 

s2  u2 


where 


is  the  total  mean  square  slope  of  the  rough  surface. 
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b .  Exponential  surface  height 
joint  probability  density 


For  the  sake  of  comparison  with  the  preceding  model,  we  select 
an  exponential  probability  density  model  of  the  following  form: 

42-2^,p(r)+t.'2'1 

(51)  pEtt.C)a— 177  r;  ni 

2ir  hz[  1  -p  2  (r)]2 


I  h  [  1-  p2(r)] 


Equation  (46)  for  J  is  solved  using  the  same  restrictions  as  with  the 
preceding  model. 


(51) 


J  = 


12 

s2u2 


\T6  t_ 
s  u 


The  optics  techniques  show  that  for  very  rough  sxirfaces,  the 
only  parameter  reflecting  the  degree  of  roughness  in  the  scattering 
cross  section  is  s,  the  rms  roughness  slope.  The  rms  roughness 
height,  h,  does  not  appear  explicitly,  and  as  such  is  not  easily 
determined  from  observations  of  average  scattered  power  from  such 
a  surface . 


2 .  Polarization  dependence 

The  elements  (3^5  will  now  be  given  for  several  polarization 
states . 


a.  Vertical  and  horizontal  states 


For  the  vertical  and  horizontal  polarization  states,  the  scattering 
matrix  elements  0vv>  Phh>  Pvh<  and  Phv  maY  be  found  upon  simplifi¬ 
cation  of  Semenov's  results,  or  from  either  of  the  other  two  approaches 

to  the  problem.  They  are: 


(53a) 


P 


vv 


sin  8j  sin  6g  sin2  4>s  Rx(i)  +  a2  a3  Rt|  (Q 

aja* 


(53b) 


Phv  =  '  sin 


az  sin  0S  (i)  -  a^  sin  9i  giijO 

al  a4 
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(53c) 


pvh  =  -  sin  4>s 


aa  sin  6S  Rx(t)  -  ag  sin  Oj  R„(i.) 

ai  a* 


(53d)  Phh  = 


sin  0^  sin  08  8in2  <))s  R|((t)  +  a2  aj  Rx(i) 


ai  84 


where  RM  (rj),  Rx(  p)  and  t  are  defined  previously  in  Eqs .  (35)  and  (36) 
and  ai>  slz,  a$,  a*  are  given  below . 


(54a) 

(54b) 

(54c) 

(54d) 


l  =  1  +  sin  bi  sin  0S  cos  4>s  -  cos  ®i  cos  0S 

a.z  -  cos  sin  08  +  sin  cos  0S  cos  4>s  t  , 

aj  =  sin  ©i  cos  08  +  cos  Q{  sin  0S  cos  <)>8, 

i 

?>4  =  COS  0|  +  COS  08 


From  the  above  expressions  and  the  two  probability  density 
models  previously  given,  one  can  substitute  the  P's  and  J  into  Eq.  (45) 
to  obtain  the  average  scattering  cross  ruction  per  unit  area. 


b.  Circular  states  and  arbitrary 
polarized  linear  Btates 

The  {3*8  given  above  for  the  vertical  and  horizontal  states,  being 
directly  proportional  to  the  scattering  matrix  elements,  can  be  sub¬ 
stituted  into  Eqs.  (27)  of  the  preceding  part  to  give  the  elements  PlR> 
Prr,  Prl,  and  PlL  *or  the  circular  polarization  states.  Equations  (29) 
and  (30)  of  that  section  can  be  used  to  give  j3r^r|g  for  linear  polarization 
states  arbitrarily  directed. 


3  .  Results  for  backscattering 

Backs catte ring  is  the  most  important  and  commonly  occurring 
situation  for  the  radar  engineer.  As  such,  the  above  results  will  be 
specialized  to  this  case,  and  curves  for  the  two  probability  density 
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models  as  a  function  of  incidence  angle  will  be  presented.  The  expres¬ 
sions  for  the  p's  given  in  Eqs .  (53)  are  allowed  to  approach  the  back- 
scattering  limit  (i.e.,  9S  -*  0i,  4>s  -*  it  )  to  yield 

(55a)  Pvv=  Phh  ~  -  sec  Qi  R(i(0)  **  sec  RJO)  , 


(55b)  Pvh  =  Phv  ""*■  ®  • 

,  * 

The  circular  polarization  backscattering  matrix  elements  found 
from  substituting  these  expressions  into  Eqs.  (27)  are 

(55c)  0LR  =  0RL  s  Pw  =  Phh  sec6iR(|(0)  , 

(55d)  PLL  =  pRR  —  0  . 

The  backscattering  matrix  element  between  arbitrary  linear 
states  as  defined  and  given  in  Eq.  (28)  becomes 

(55e)  pr^qg  =  Pw  cos^s  -  m)  a  -  sec  0i  Rn(0)  cos^g  -  r\{)  . 

The  above  equations  demonstrate  the  claim  of  Hagfors  and  other 
investigators  who  note  that  optics  techniques  applied  to  very  rough 
surfaces  models  predict  no  depolarization  for  backscattering.  The 
second  and  third  optics  models  discussed  previously  clearly  show  that 
backscattering  can  come  only  from  areas  or  facets  on  the  rough  surface 
which  are  oriented  normal  to  the  incidence  direction.  This  also 
accounts  for  ^he  appearance  of  the  Fresnel  reflection  coefficients  for 
normal  incidence  in  the  above  equations . 

Let  us  now  give  the  results  for  average  backscattering  cross 
sections-  per  unit  area  for  the  two  probability  density  models  proposed 
above,  remembering  that  no  depolarization  is  predicted. 
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a.  Gaussian  surface  height  probability 
function  model 


(56) 


vv 


=  cr 


o  _ 


hh 


=  O' 


RL 


=  cr 


LR 


=  a 


Hi’li 


_  sec4^ 
s2 


v°>l 


tan2 

s2 


b.  Exponential  surface  height 
probability  function  model 


(57) 


vv 


<r 


o 


3  sec4^ 
s2 


tan 


The  first  of  the  above  equations  has  been  given  previously  by 
Hagfors  (1964)  while  the  second  is  presented  here  for  comparison. 
Curves  illustrating  the  dependence  of 


|r,0)I2  i 

for  the  above  two  models  upon  angle  of  incidence,  ^i,  for  various  values 
of  rms  roughness  slope,  s,  are  presented  in  Fig.  6. 


IV.  COMPOSITE  SURFACES 
»  ^  A  Explanation 

I 

*  Natural  forces  seldom  create  a  "very  rough"  surface  satisfying 

f.  the  restrictions  of  the  preceding  part.  Rough  surfaces  arising  in  nature 
having  roughness  whose  height  is  large  compared  to  wavelength  almost 
always  also  possess  smaller  scale  roughnesses,  such  that  the  surface 
is  not  locally  smooth.  As  examples,  consider  an  area  of  mountainous 
terrain.  The  mountains  and  valleys  alone  comprise  a  roughness  large 
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Angle  From  Virtual ,  dpgrm  Angle  From  vertical ,  degrees 

Fig.  6.  Average  incoherent  backscattering  cross  section  per 
unit  area  for  very  rough  surface  statistical  model 
vs.  incidence  angle  for  any  set  of  direct  polarization 
states  for  two  surface  height  JPDF  (Joint  Probability 
Density  Function)  models .  Curves  are  normalized  to 
the  Fresnel  reflection  coefficients  for  normal  incidence. 


compared  to  wavelength  for  frequencies  of  VHF  or  higher.  If  no  other 
smaller  roughness  were  present,  the  analysis  of  the  preceding  section 
might  be  applied.  However,  the  trees,  rocks,  and  grassy  fields  also 
present  constitute  a  roughness  whose  scale  compared  to  wavelength 
may  be  small.  Hence,  the  tangent  plane  restriction  and  optics  approxi¬ 
mation  of  the  preceding  part  cannot  be  used  without  closer  examination. 
The  type  of  surface  described  above,  characterized  by  both  large  and 
small  scale  roughnesses,  will  be  called  a  "composite"  surface. 

Let  us  look  again  at  the  chief  differences  between  the  average 
incoherent  backscattered  power  from  the  slightly  rough  surface  and 
the  very  rough  surface  of  the  preceding  parts.  The  slightly  rough 
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surface,  as  shown  in  Fig.  2,  produces  a  small  backscattered  field  near 
normal  incidence  where  0^  *  0,  but  this  field  may  still  be  significant 
near  grazing.  As  discussed  previously  the  scattered  power  near  graz¬ 
ing  is  directly  related  to  the  amount  of  surface  roughness  spectral 
strength  at  higher  spatial  frequencies;  in  particular,  at  and  near 
frequency  2k0  radians/meter.  On  the  other  hand,  backscattering 
from  very  rough  surfaces,  as  seen  in  Fig.  6,  is  much  larger  near 
normal  incidence  (0i  3  0)  but  falls  off  rapidly  away  from  normal. 

The  suggestion  to  be  put  forth  here  is  that  both  the  large  scale 
and  the  small  scale  roughness  contribute  significantly  to  the  scattering 
from  natural,  composite  surfaces.  Near  the  specular  direction  (i.e., 

0S  a  ®i»  ^s  s  0)  scattering  of  the  optics  type,  predicted  by  the  very 
rough  surface  theories,  predominates.  At  scattering  angles  consider¬ 
ably  removed  from  this  direction,  however  the  smaller  but  persistent 
scattered  power  is  due  to  the  slight  roughness  present  on  top  of  the 
larger  roughnesses. 

In  fact,  to  a  first  approximation,  one  may  merely  add  the 
average  incoherent  scattering  cross  sections  from  the  very  rough 
surface  model  to  that  of  the  slightly  rough  surface  model  to  obtain 
the  total  composite  rough  surface  scattering  cross  section.  An 
heuristic  and  physically  intuitive  proof  of  this  statement  will  be  given 
here,  but  for  a  more  exact-albeit  mathematically  involved-proof,  one 
can  consult  two  recent  Soviet  articles  of  Semenov  (1966)  and  Fuks  (1966). 

Consider  first  the  very  rough  surface  model  of  the  preceding 
part.  The  surface  to  which  it  applies  is  smooth  over  regions  of  the 
order  of  a  wavelength.  The  scattered  power  is  incoherent  and  con¬ 
centrated  near  the  specular  direction  because  the  surface  slopes  are 
generally  small  and  the  probability  is  great  that  the  local  normal  of 
any  surface  point  is  within  a  few  degrees  of  the  vertical.  Now  assume 
that  a  slight  roughness  "grows"  on  top  of  this  very  rough  but  locally 
smooth  surface.  This  slight  roughness  does  not  greatly  diminish  the 
large  specular  scattering  of  the  very  rough  surface,  just  as  it  did  not 
diminish  greatly  the  coherent  scatter  of  the  smooth  plane  over  which 
it  was  imposed  in  Part  2.  However,  in  directions  considerably  away 
from  the  specular  (e.g.,  near  grazing  for  backscatter),  where  the 
very  rough  surface  model  predicts  very  small  scattered  power,  the 
slight  roughness  will  produce  an  observable  addition  to  the  scattered 
power.  Since  the  "very  rough"  portion  of  the  surface  has  small  sur¬ 
face  slopes,  the  slight  roughness  can  be  approximated  to  a  first  order 
as  being  distributed  over  a  smooth  plane  in  computing  the  local  angle 
of  incidence.  As  seen  in  Figs.  2  and  3,  small  changes  (say  15*  or 
less)  in  the  local  angle  of  incidence,  0*,  (due  to  local  variations  in  the 
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underlying  very  rough  surface  slope)  do  not  produce  a  marked  change 
in  backscattering  cross  section,  so  long  as  the  spectral  strength  of 
the  superimposed  slight  roughness  is  fairly  constant  out  to  approxi¬ 
mately  i  kD  radians/meter.  Hence  the  two  average  incoherent  scat¬ 
tering  cross  sections  may  be  simply  added  together,  assuming  the 
same  incidence  and  scattering  angles  (B,,  9s,4>s)  f°r  both. 

Actually,  two  averaging  processes  must  take  place.  The  slight 
roughness  height  variable  and  the  large  roughness  height  variable  are 
averaged  separately.  This  necessitates  the  additional  restriction  that 
the  two  different  scales  of  roughness  are  statistically  independent. 

Such  is  generally  true  when  the  two  roughness  scales  are  formed  by 
different  natural  processes.  For  example,  mountains  and  valleys 
are  usually  formed  by  an  entirely  different  process  than  is  the 
vegetation  cover;  hence  the  assumption  of  statistical  independence 
between  the  distribution  of  these  natural  roughnesses,  although  not 
perfect,  generally  seems  reasonable.  Since  the  slight  roughness  and 
large  roughness  are  statistically  independent  and  each  alone  produces 
strictly  an  incoherent  scattered  power,  we  can  perform  the  averaging 
of  each  separately  and  merely  add  the  powers  from  each. 

The  above  explanation  is  meant  to  be  intuitive  rather  than 
rigorous-  It  is  meant  to  provide  physical  insight  rather  than  mathe¬ 
matical  confusion.  Mathematical  proof  is  satisfying,  however,  and 
can  be  found  in  the  above  mentioned  Soviet  references.  Generally 
the  necessary  restrictions  are  those  given  in  the  preceding  two  parts, 
along  with  the  one  mentioned  in  the  preceding  paragraphs.  Neglected 
in  this  explanation  is  any  possible  roughness  scale  which  is  of  inter¬ 
mediate  size  such  that  it  does  not  fall  entirely  either  in  the  very 
rough  or  the  slightly  rough  category.  Since  most  natural  surfaces 
have  a  roi^ghness  component  whose  scale  is  of  this  order,  the  model 
here  musij  be  expected  to  give  results  which  deviate  somewhat  from 
measured'  values .  This  deviation  is  expected  for  scattering  geometries 
where  the  scattering  cross  section  for  the  slight  roughness  alone  is 
approximately  equal  to  that  for  the  large  scale  roughness.  In  general 
the  predicted  result  will  be  too  small,  because  any  intermediate  scale 
roughness  will  tend  to  add  another  scattered  field  component  which 
will  be  significant  in  this  region.  Since  no  known  theory  is  valid  for 
intermediate  scale  roughness,  this  component  cannot  at  present  be 
taken  into  account  quantitatively. 

Western  writers  have  attempted  to  treat  composite  rough  sur¬ 
faces .  Beckmann  (1965)  and  Hayre  and  Kaufman  (1965)  both  treat  this 
problem  by  a  physical  optics  model.  Their  analysis  does  not  truly 
apply  to  a  slight  roughness  scale  since  it  is  based  on  an  optics  formulation. 
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This  is  evident  from  the  results,  which  are  functionally  the  same  as 
for  the  very  rough  surface  with  a  single  large  scale  of  roughness.  The 
only  difference  is  that  the  mean  square  slope,  s2,  appearing  in  their 
results  for  the  composite  surface  is  written  in  terms  of  the  sum  of  the 
mean  square  slopes  of  each  of  the  roughness  scales  present.  Their 
results  are  hence  valid  for  all  roughness  scales  down  to  those  which 
can  no  longer  be  treated  by  their  optics  formulation  (i.e.,  which  no 
longer  satisfy  the  restrictions  of  Part  3).  They  cannot  apply  to  sur¬ 
faces  which  include  a  true  slight  roughness  component,  however. 

Radar  experimentalists  have  in  the  past  attempted  to  break  the 
measured  direct  polarized  backscatter  cross  section  vs.  incidence 
angle  into  two  regions.  The  region  near  normal  was  called  the 
specular  region,  and  it  was  generally  agreed  to  be  predicted  by  the 
optics  theories  and  results  of  Part  3.  The  "tail"  of  the  echo,  or  that 
region  near  grazing  was  termed  the  "diffuse"  region,  and  there  has 
been  no  general  agreement  to  date  on  a  satisfactory  explanation  for 
its  existence.  Various  investigators  have  attempted  to  fit  empirical 
curves  to  this  component;  among  the  postulated  are  cos  0i,  cos3/2  0(, 
and  still  others.  Hagfors  (1967)  has  suggested  a  model  consisting  of 
randomly  oriented  dipoles  as  accounting  for  this  component.  This 
is  one  of  a  large  class  of  geometrical  models  with  distributions  of 
individual  scatterers  which  lie  outside  the  class  of  continuous  sur¬ 
faces  we  have  been  considering.  (See  Peake,  1967). 

We  wish  to  suggest  here  that  this  "diffuse"  component  is  pro¬ 
duced  by  the  presence  of  slight  and  intermediate  roughness  scales 
and  is  predicted  to  a  first  order  by  the  models  of  Part  2.  We  offer 
predicted  curves  based  on  the  composite  model  presented  here  and 
compare  them  with  measured  results  for  backscatter. 

The  explanation  for  the  depolarized  backscattered  component 
from  composite  surfaces  is  more  complex,  however.  The  optics 
theories  predict  no  return.  To  the  first  order,  the  slightly  rough 
model  studied  here  also  indicates  no  depolarized  return  for  the  verti¬ 
cal  and  horizontal  states;  depolarization  for  the  circular  states  does 
appear,  but  is  absent  at  normal  incidence.  The  higher  order  terms, 
for  example  the  second  term  as  calculated  by  Valenzuela,  do  show 
depolarization,  even  for  the  vertical  and  horizontal  states.  As  the 
roughness  scale  increases  with  respect  to  wavelength  from  slight  to 
intermediate,  these  higher  order  terms  containing  depolarization 
effects  become  non- negligible .  Since  the  perturbation  theory  breaks 
down  as  the  smallness  parameters  approach  unity,  the  contribution 
of  the  intermediate  scale  of  roughness  cannot  be  calculated.  It  is  clear 
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from  the  trend  of  the  higher  order  terms,  however,  that  they  will  con¬ 
tribute  significantly  to  both  the  polarized  and  depolarized  components. 

It  is  our  conclusions,  therefore,  that  depolarization  for  backscatter  is 
produced  by  two  mechanisms:  (a)  the  slight  and  intermediate  roughness 
scales  just  discussed,  and  (b)  multiple  reflections  from  suitably 
oriented  surface  elements  of  the  large  scale  roughness  for  surfaces 
where  the  reflection  coefficients  are  polarization  dependent. 

To  illustrate  the  predictions  of  the  composite  surface  model, 

Figs.  7  provide  curves  for  horizontal  and  vertical  polarization  states 
for  a  variety  of  dielectric  constants  and  model  parameters.  They  are 
made  by  simply  adding  the  backscattering  cross  sections  for  the  slightly 
rough  model  to  those  of  the  very  rough  model.  More  curves  of  this  type 
can  be  found  in  Radar  Cross  Section  Handbook  (1968).  Figure  8  shows 
the  measured  cross  section  of  the  sea  surface  on  a  calm  day  made 
for  the  vertical  and  horizontal  states;  these  are  also  taken  from  the 
above  reference.  The  similarity  is  convincing.  Note  that  none  of 
these  curves  represents  a  depolarized  component;  the  received 
polarization  state  is  in  each  case  the  same  as  the  transmitted. 

For  the  circular  polarization  states,  one  can  construct  the 
same  type  of  curves,  which  results  are  sort  of  a  superposition  of  the 
curves  of  Fig.  3  with  those  of  Fig.  6.  They  are  not  shown  here  but 
may  be  found  in  the  preceding  reference.  It  should  be  noted,  however, 
that  the  curve  representing  the  depolarized  component  is  not  actually 
a  composite,  because  the  very  rough  surface  contribution  to  back- 
scattering  in  this  case  is  zero:  only  the  slight  roughness  produces  de¬ 
polarization.  One  set  of  curves  for  circular  polarization  is  compared 
to  a  measured  set  from  the  lunar  surface  at  \Q  =  68  cm  in  Fig.  9. 
Measured  results  were  taken  from  the  work  of  Evans  and  Pettingill 
(1963),  and  were  normalized  with  respect  to  the  measured  cross  section 
of  the  entire  moon  at  this  frequency.  The  set  of  model  curves  was 
selected  on  the  basis  of  best  fit.  From  the  model  for  this  best  fit,  the 
following  parameters  for  the  lunar  surface  are  indicated:  «r  =  dielectric 
constant  =  2.9,  tan-18i  *  large  scale  roughness  slope  =  12°,  h2  =  rms 
slight  roughness  height  =  2  cm,  tan_1s2  =  slight  roughness  slope  =  27°  • 
The  comparison  of  the  measured  and  predicted  depolarized  components 
shows  lack  of  agreement  near  normal  incidence .  This  is  expected  for 
the  following  reasons:  (i)  Hagfors  (1967)  reports  that  the  ability  to  iso¬ 
late  one  measured  circular  state  from  the  other  was  probably  not 
greater  than  18  dB.  Hence  near  normal  where  the  polarized  component 
is  very  strong,  most  of  the  "measured"  depolarized  component  may 
actually  be  the  polarized  component,  (ii)  The  intermediate  scale  rough¬ 
ness,  which  cannot  be  accounted  for  by  the  theory,  most  likely  produces 
depolarization.  As  mentioned  previously,  this  effect  is  likely  to  appear 
further  away  from  grazing,  i.e.,  near  normal  incidence. 
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(a) 

Average  backscattering  cress  section  per  unit  area  for 
composite  rough  surface  model  vs.  incidence  angle  for 
various  large  scale  roughness  parameters,  sj.  and 
small  scale  roughness  parameters,  kcf  and  kDh  (or  s2) 
Linear  polarization  states,  (a)  Gaussian  large  scale 
roughness  height  probability  density  model. 
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Fig.  7b.  Exponential  large  scale  roughness 
height  probability  density  model. 
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Fig.  8.  Measured  average  back- 
scattering  cross  section 
per  unit  area  from  re¬ 
latively  calm  sea  surface 
at  X-band  vs.  incidence 
angle  for  vertical  and 
horizontal  polarization 
states . 


Angit  From  Vertical  .degrees 
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V.  ROUGH  SPHERICAL  SURFACE  - 

BACKSCAT^ERING  CROSS  SECTION 


A.  Introduction 


The  past  sections  have  derived  the  average  incoherent  scattering 
cross  section  per  unit  area  for  a  rough  planar  surface,  i.e.,  a  surface 
whose  mean  value  is  a  plane.  In  this  section  we  will  extend  these 
results  to  a  rough  spherical  surface.  Such  a  situation  is  representa¬ 
tive,  for  example,  of  a  planetary  surface,  or  of  many  spherical 
satellites  of  the  Echo  I  and  II  variety  whose  surfaces  are  somewhat 
rough.  In  defining  the  average  backscattering  cross  section  of  a  rough 
spherical  surface,  we  understand  that  the  entire  sphere  is  immersed 
simultaneously  in  the  illuminating  field.  This  may  not  be  the  case,  for 
example,  in  many  of  the  planetary  radar  studies  where  a  short  pulse 
sweeps  past  the  planet,  illuminating  only  an  annular  area  at  a  time. 
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Fig.  9.  Measured  average  back- 
scattering  cross  section 
per  unit  area  from  lunar 
surface  at  UHF  (X.  =  68  cm) 
vs  .  incidence  angle  with 
circular  polarization  states. 
Dashed  curves  represent 
composite  rough  surface 
model,  presented  for  com¬ 
parison.  Upper  curves 
represent  opposite  sense 
polarization  states,  lower 
curves  -  same  sense. 


If 


The  results  of  the  preceding  sec.  ions  can  be  used  to  obtain  the 
incoherent  backscattering  cross  sections  for  both  very  rough  and  slightly 
rough  spheres.  Furthermore,  for  the  slightly  rough  sphere,  there  will 
be  a  significant  coherent  component  due  predominantly  to  reflection  from 
the  front  cap,  and  this  will  also  be  estimated. 

It  is  assumed  that  the  radius  of  the  sphere  Ar,  is  much  larger 
than  both  the  wavelength,  X,  and  roughness  height  correlation  length, 
i  .  It  is  further  assumed  that  the  sphere  is  large  enough  and  its 
material  lossy  enough  that  waves  entering  the  sphere  do  not  further 
contribute  to  the  scattering.  In  addition,  all  of  the  restrictions  of  the 
preceding  sections  apply.  j 
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B  •  Coherent  Backscattering  Cross  Section 


It  was  mentioned  in  Part  I  that  the  coherent  scatter  from  any 
slightly  rough  surface  is  highly  dependent  upon  the  shape  of  the  illumi¬ 
nated  surface  in  the  absence  of  roughness.  The  presence  of  roughness 
is  taken  into  account  by  modifying  the  Fresnel  reflection  coefficients  at 
the  surface,  as  in  Eq.  (1).  The  reflected  complex  fields  in  this  case, 
being  coherent  and  having  a  non- zero  average  value,  can  be  treated 
just  as  they  would  be  from  a  smooth  surface;  stationary  phase  analysis 
predicts  that  nearly  the  entire  coherent  backscattered  field  in  this  case 
comes  from  the  front  specular  cap.  Hence  the  average  coherent  back- 
scatter  cross  section  from  a  slightly  rough  spherical  surface  becomes 

2  2 

(58)  <r  =  ir  |  R(0)  | z  e'4k°  h 
where 

n/TT- 

(59)  R(0)  =  — 

v«r  + 


is  the  Fresnel  reflection  coefficient  for  a  smooth  planar  surface  at 
normal  incidence.  The  effect  of  surface  roughness  near  the  specular 
cap  is  accounted  for  by  the  exponential  factor  involving  the  mean  square 
roughness  height,  h2.  The  above  result,  having  come  from  the  use  of 
the  slightly  rough  surface  reflection  coefficients  of  Eq.  (1)  is  valid  for 
surfaces  having  a  Gaussian  distributed  roughness  height.  One  sees 
that  as  rms  roughness  height  increases  to  where  kDh  >  1,  the  coherent 
cross  section  decreases  rapidly. 

The  polarization  states  of  the  scattered  coherent  field  above  is 
the  same  as  it  would  be  in  the  absence  of  roughness,  i.e.,  there  is  no 
depolarization  predicted. 


C.  Average  Incoherent  Backscattering  Cross 
Section  -  Slightly  Rough  Surface 

The  average  incoherent  backscattering  cross  sections  for  a 
slightly  rough  spherical  surface  are  computed  from  the  preceding 
results  for  a  rough  planar  surface.  Use  is  made  of  the  property  of 
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incoherent  power  which  permits  one  to  merely  add  the  power  backseat  - 
tered  from  one  portion  of  the  surface  to  that  from  another  portion, 
rather  than  adding  the  complex  scattered  fields  themselves.  Recall 
thatcr^  is  the  average  incoherent  backscattered  power  per  unit  mean 
surface  area  (where  9S  =  6;  and  4>s  =  it  in  the  equations  for  cr  ^ g) .  Then 
the  backscattering  cross  section  for  a  mean  spherical  area  element 
dA  =  sin  d9i  d4>  is  the  product  of  dA  and  cr  °g .  The  total  average 
backscattering  cross  section  is  then  determined  Dy  integrating  over  the 
entire  illuminated  hemisphere,  i.e.. 


(60) 


V6 


2 

=  ar 


2tt 


I 


tt/2 


y6 


sin  9i  d6i 


where  the  integration  over  <j>  has  been  performed  since  the  integrand  is 
not  a  function  of  <j> .  The  spherical  coordinate  system  has  been  chosen 
here  so  that  the  polar  axis  coincides  with  the  direction  of  propagation, 
making  the  polar  angle  9^  also  the  angle  of  incidence  upon  the  spherical 
surface . 


For  the  aligned  and  crossed  incident  and  received  linear  polari¬ 
zation  states,  the  expressions  which  must  be  substituted  into  the  above 
equation  for  o-  are  those  given  in  Eqs.  (19)  and  (31).  The  integration 
over  <j>  corresponds  to  the  integration  over  which  was  already  per¬ 
formed  in  deriving  Eqs.  (31).  The  integration  over  9i  demanded  in  the 
above  equation  can  only  be  performed  numerically.  The  resulting  cross 
sections  normalized  as 


_ a _ 

»  A  £  h2 

are  plotted  in  Fig.  10  vs  koi  for  various  values  of  surface  dielectric 
constant.  This  is  done  for  both  the  aligned  and  crossed  linear  states 
(representing  the  polarized  and  depolarized  return)  and  for  the  Gaussian 
and  exponential  surface  height  correlation  coefficient  models.  Note  that 
as  kof  becomes  larger  the  depolarized  component  becomes  smaller. 

This  is  expected,  since  larger  koi  corresponds  to  a  smoother  surface. 
Note  that  the  terminology  vertical  and  horizontal  states  have  no  signifi¬ 
cance  when  referred  to  a  sphere. 

For  the  circular  polarization  states,  the  appropriate  curves  are 
found  in  Radar  Cross  Section  Handbook  (1968).  They  are  nearly  identical 
with  those  for  the  linear  states  in  the  cases  of  the  polarized  components  { 
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(a) 


Fig.  10a.  Average  incoherent  backscattering  cross  section  for 
slightly  rough  spherical  surface  of  mean  radians  Ar 
as  a  function  of  koi.  Cross  section  is  normalized  to 
tt  Ar2k02h2.  Linear  polarization  states,  upper  curves 
represent  same  sense  received  as  transmitted,  lower 
curves  -  perpendicular  sense  received,  (a)  Gaussian 
height  correlation  coefficient  model. 
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I  Fig.  10b.  Exponential  height  correlation  coefficient  model. 

I 
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the  depolarized  component  for  the  circular  states  is  exactly  twice  that 
for  the  depolarized  linear  states,  as  seen  from  Eqs .  (31b)  and  (28b). 


D.  Average  Incoherent  Backscattering  Cross 
Section  -  Very  Rough  Surface 

The  procedure  to  be  followed  to  determine  the  backscattering 
cross  section  of  a  sphterical  surface  when  the  irregularity  falls  into  the 
very  rough  category  is  the  same  as  that  of  the  preceding  section.  The 
integral  of  Eq.  (60)  is  solved  (again  numerically)  using  or  °  predicted 
in  Eqss.  (56)  and  (57)  of  Part  III.  Recall  that  the  backscattered  power 
predicted  in  this  case  is  not  depolarized,  and  the  polarized  component 
is  the  same  for  the  aligned  linear  states  as  for  the  circular  states . 
Hence  only  one  set  of  curves  need  be  plotted  here. 

The  average  backscattering  cross  sections  are  shown  in  Fig.  11, 
normalized  as  °‘/(ttAr2  |  R (0)  1 2  ) 


tan'i.aeqreM  (s»,ms  rouQfwtt* *iope * y) 

Fig.  11.  Average  incoherent  normalized  backscattering 
cross  section  for  a  very  rough  spherical  sur¬ 
face  of  mean  radius  Aj^  as  a  function  of  tan"ls, 
where  s  is  the  rms  roughness  slope  cross 
section  is  normalized  to  tt  Ar2|  R(0)  |  2,  where 
R(0)  is  the  Fresnel  reflection  coefficient  for 
normal  incidence. 
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and  plotted  as  a  function  of  the  arctangent  of  the  rms  roughness  slope, 
s  (i.e.,  s  =  (2h/i)).  This  function  is  sometimes  referred  to  as  the 
roughness  gain  of' the  sphere.  The  two  curves  represent  the  two 
probability  density  function  models  discussed  in  Part  III.  Noteworthy, 
however,  is  the  fact  that  as  roughness  increases,  the  backscattering 
cross  sections  become  larger  than  that  for  a  smooth  sphere. 

It  is  seen  that  as  the  surface  becomes  smooth  (i.e.,  s  — *•  0),  the 
backscattering  cross  sections  do  indeed  approach  those  for  a  smooth 
sphere  i.e.,  a  =  ir  |R(0)|2  .  The  figures  are  of  questionable . 
validity  for  tan-1  s  >  45°  because  of  the  restrictions  under  which 
they  were  derived.  For  typical  surface  slopes  much  beyond  45*, 
shadowing  and  multiple  scattering  are  no  longer  negligible.  In  the 
roughness  slope  region  less  than  45®  the  increase  in  scattering  cross 
section  from  the  smooth  sphere  limit  is  not  appreciable.  Hence  an 
important  conclusion  for  planetary  surface  backscatter  caji  be  drawn: 

Large  scale  surface  roughness  for  all  practical  purposes 
can  be  neglected  and  the  total  average  backs catteyii^g  cross 
section  can  be  taken  to  be  that  for  a  smooth  sphere,' 

o-  =  itA^  I  R(0)  | 2  .  /  1 

!  i 

Consequently,  knowing  the  approximate  radius  of  the  body,  Ar,  an 
estimate  of  the  surface  dielectric  constant  can  be  made  from  a 
measurement  of  cr  /  tt  Aj^  ;  for  the  moon  for  the  range  5  cm  '<  X.  <100  cm 
one  obtains  approximately  |r(0)|2  =  0.07.  Such  a  value  of  J  R( 0)  [ 2 
implies  an  effective  dielectric  constant  of  cr  =  2.9.  \ 

\  t 

E .  Average  Incoherent  Backscattering  Cross 
Section  -  Composite  Surface 

As  mentioned  previously,  most  natural  surfaces  consists  of 
roughness  scales  both  larger  and  smaller  than  wavelength.  As  shown 
in  Part  IV,  one  can  add  the  incoherent  backscattering  cross  sections  of 
the  preceding  sections  (presented  in  Figs.  10  and  11)  to  obtain  the 
effective  total  backscattering  cross  section  of  the  composite  surface. 
However,  the  backscattering  cross  section  for  the  slight  roughness, 
being  proportional  to  k02h2,  is  small,  since  this  parameter  is  re¬ 
stricted  to  small  values.  Hence  only  the  large  scale  roughness  contributes 
significantly  to  the  backscattered  power  for  the  polarized  received  com¬ 
ponent.  On  the  other  hand,  since  only  the  slight  roughness  scale  is 
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responsible  for  the  depolarized  component,  according  to  these  models, 
it  alone  may  be  the  best  available  estimate  for  this  component. 


VI.  SUMMARY 

We  have  attempted  to  set  forth  what  we  believe  to  be  the  soundest 
analyses  and  solutions  for  rough  surface  scattering  that  are  presently 
available.  In  addition,  we  have  enumerated  the  restrictions  under  which 
the  results  are  valid.  Finally,  we  have  interpreted  physically  the  scat¬ 
tering  mechanism  behind  the  mathematics. 

.  A  slightly  rough  surface  produces  strong  coherent  reflection 
near  the  specular  direction  and  a  weaker  incoherent  scatter  in  other 
directions.  An  estimate  of  the  coherent  reflection  can  be  computed 
just  as  it  is  fov  a  smooth  surface,  making  use  of  the  reflection  coef¬ 
ficients  of  Eq.  (1);  the  resulting  scattering  pattern  is  strongly  de¬ 
pendent  upon  the  shape  of  the  illuminated  area  and  exhibits  a  lobe 
structure.  The  incoherent  scatter  was  dealt  with  in  more  detail  by 
a  perturbation  technique  rather  than  the  less  valid  tangent  plane  approxi¬ 
mation.  It  was  shown  that  incoherent  scattered  power  from  a  slightly 
rough  surface  of  either  homogeneous  material  or  a  perfectly  reflecting 
interface  is  directly  proportional  to  the  roughness  spectral  densities. 

The  highest  roughness  spectral  components  which  can  affect  the  process 
are  those  near2k0  radians/meters.  These  are  responsible  for  back- 
scatter  near  grazing  incidence.  Lower  roughness  spectral  frequencies 
produce  scattering  closer  to  the  specular  direction.  For  the  important 
case  of  backscatter,  significant  differences  exist  between  the  cross 
sections  near  grazing  for  the  vertical  and  horizontal  states;  vertical 
transmit  and  receive  antennas  produce  considerably  more  return  than 
horizontal,  a  fact  which  is  confirmed  experimentally.  The  solutions, 
including  only  the  lowest  order  perturbation  terms,  exhibit  no  de¬ 
polarization  for  the  vertical  and  horizontal  states.  From  the  solutions 
for  vertical  and  horizontal  polarizations,  results  for  the  circular  states 
and  arbitrarily  oriented  linear  transmit  and  receive  states  are  given. 

The  only  presently  satisfactory  methods  for  treating  very  rough 
surfaces  are  optics  techniques;  basic  to  all  of  these  is  the  tangent 
plane  restriction.  A  physical  interpretation  of  the  scattering  process 
can  be  based  upon  any  of  the  well  known  optics  principles.  Three 
techniques  (i.e.,  physical  optics,  ray  optics,  and  the  specular  point  or 
stationary  phase  principle)  all  give  identically  the  same  result.  Scatter¬ 
ing  in  a  given  direction  takes  place  as  a  result  of  surface  facets  so  oriented 
that  they  specularly  reflect.  As  would  be  expected  from  such  a  theory, 
backscattered  power  is  not  depolarized.  Two  probability  density  function 
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models  are  chosen  for  the  surface  height;  for  both,  the  curves  of  b>».ck- 
scattered  power  show  strong  return  near  normal  incidence  (the  specular 
direction)  which  falls  off  rapidly  as  grazing  is  approached.  The  only 
surface  height  correlation  coefficient  which  can  be  employed  with  these 
techniques  is  one  which  is  parabolic  in  behavior  for  small  surface 
separations . 

Natural  surfaces  which  are  very  rough  in  height  compared  to 
wavelength  almost  always  have  smaller  scale  roughnesses  superimposed; 
these  are  called  composite  surfaces.  As  one  might  expect,  the  slightly 
rough  and  very  rough  surface  scattering  theories  together  account  in 
part  for  the  return.  In  the  case  of  backscatter  near  vertical,  for 
instance,  the  return  is  dominated  by  the  very  rough  surface  predicted 
scatter.  Near  grazing,  however,  the  smaller  roughness  scales  usually 
account  for  the  return.  Curves  for  backscattering  cross  section  are 
shown  for  these  composite  models  at  various  set  of  parameters  for  the 
vertical  and  horizontal  polarization  states.  These  are  compared  with 
measured  data  from  the  lunar  and  sea  surface. 

Backscatter  from  roughened  spherical  surfaces,  which  has 
application  to  planetary  probing  and  passive  satellite  communitations, 
is  treated  in  Part  V.  The  coherent  and  incoherent  average  backscat¬ 
tering  cross  sections  are  obtained,  the  latter  in  graphical  form  resulting 
from  a  numerical  integration.  Very  rough  spherical  surfaces  are  also 
treated,  and  it  is  shown  that  the  roughness  in  this  case  has  little  effect 
on  the  backscatter  relative  to  that  of  a  smooth  sphere;  in  reality,  the 
roughness  slightly  increases  the  cross  section. 
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